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Abstract 

The Regge behaviour of the scattering amphtudes in perturbative QCD is governed in 
the generahzed leading logarithmic approximation by the contribution of the color-singlet 
compound states of Reggeized gluons. The interaction between Reggeons is described by 
the effective hamiltonian, which in the multi-color limit turns out to be identical to the 
hamiltonian of the completely integrable one-dimensional XXX Heisenberg magnet of non- 
compact spin s = 0. The spectrum of the color singlet Reggeon compound states, - per- 
turbative Pomerons and Odderons, is expressed by means of the Bethe Ansatz in terms of 
the fundamental Q— function, which satisfies the Baxter equation for the XXX Heisenberg 
magnet. The exact solution of the Baxter equation is known only in the simplest case of 
the compound state of two Reggeons, the BFKL Pomeron. For higher Reggeon states the 
method is developed which allows to find its general solution as an asymptotic series in 
powers of the inverse conformal weight of the Reggeon states. The quantization conditions 
for the conserved charges for interacting Reggeons are established and an agreement with 
the results of numerical solutions is observed. The asymptotic approximation of the energy 
of the Reggeon states is defined based on the properties of the asymptotic series, and the 
intercept of the three-Reggeon states, perturbative Odderon, is estimated. 
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1. Introduction 



Recently, a lot of attention has been renewed to the old problem of understanding QCD Pomeron. 
The interest was partially inspired by new exciting experimental results, which confirmed a 
growth of the structure function of deep inelastic scattering, F2{x,Q'^), at small values of the 
Bjorken variable x and increasing with energy of the total hadronic cross-sections, atot{s)- Both 
phenomena allow us to test QCD in the extreme limit of high energies and fixed transferred 
momenta, the limit in which the famous Rcgge model emerges [1]. Being one of the most popular 
subjects in 70's, the Regge theory interprets the increasing of physical cross-sections with the 
energy by introducing the concept of the Regge poles. The Regge poles can be classified into 
different families according to their quantum numbers and their contribution to the amplitude of 
elastic scattering of hadrons hi and /i2 in the Regge limit, s 3> —t, has the form A/jj/j2(s, t) ~ s"^*^ , 
with a = a{t) being universal Regge pole trajectory independent on particular choice of the 
scattered hadrons [1]. 

Among all possible families of the Regge poles there is a special one, with the quantum 
numbers of vacuum, the so-called Pomeron. The Pomeron provides the leading contribution to 
the higher energy behaviour of the total hadronic cross-section [1] 

<hM = s-'lmA^,,,{s,t = 0) ~ 

with ap(0) being the Pomeron intercept. Replacing one of incoming hadrons by virtual photon 
with invariant mass —Q^ we may find a similar expression for the asymptotics of the structure 
function of deep inelastic scattering YiQ) + h —>■ X in the Regge limit, x = Q'^/s <C 1 and 

= fixed. There is however an important difference between these two processes from point 
of view of underlying QCD dynamics. The Regge behavior of the total hadronic cross-sections 
is essentially nonpcrturbativc phenomenon and it is governed by the "soft" QCD Pomeron [2] 
with the intercept 0:^(0) ~ 1.08. At the same time, the deep inelastic scattering can be analyzed 
for large enough within the framework of perturbative QCD by using the operator product 
expansion, or factorization, and the small— x asymptotics of the structure function is controlled 
by the "hard" Pomeron [2], F2{x,Q'^) ~ x^-"ip(°\ with larger intercept q;i>(0) ^ 1.38. 

Although the Regge model is extremely successful in describing a wealth of different data 
[1, 2] it is still unclear whether it can be justified from the first principles of QCD. Since we 
do not have a complete understanding of nonpcrturbativc QCD we have to restrict ourselves 
to processes involving "hard" Pomeron. As an example of such processes, one may consider 
the scattering of two onium states each consisting of a heavy quark-antiquark pair [3, 4]. The 
first attempts to describe the "hard" Pomeron in QCD revealed the remarkable property of 
gluon reggeization [3, 5]. It was found from the analysis of the Feynman diagram in the leading 
logarithmic approximation (LLA), aglns ~ 1 and as <^ 1, that interacting with each other 
gluons form a new collective excitation, Reggeon, which according to its contribution to the 
scattering amplitude can be interpreted as a Regge pole with quantum numbers of a gluon. As a 
result, QCD can be replaced in the Regge limit by an effective field theory [6] in which Reggeons 
play a role of new elementary fields while Pomerons appear as compound states of interacting 
Reggeons. The simplest example of such state built from only two Reggeons, Balitsky-Fadin- 
Kuraev-Lipatov Pomeron, was found many years ago [3] and its intercept was calculated in the 
LLA as 

a N 

«BPKL = l + ^41n2 (LI) 
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where as is the QCD coupHng constant and N is the number of quark colors. Apart from the 
BFKL Pomeron there is an infinite number of another states with vacuum quantum numbers, 
which arc built from an arbitrary number of Rcggeons n > 2. 

To identify higher n Reggeon compound states we have to go beyond the LLA in calculating 
the scattering amplitudes. One possible way of doing this has been proposed in [7, 8] and is called 
the generalized leading logarithmic approximation. In this approximation, the amplitude of the 
onium-onium scattering, A{s,t), is given in the Regge limit by the sum of diagrams describing 
the propagation of an arbitrary number n of interacting Reggeons in the channel 

oo 

A{s,t) = J2<-'Ms,t). (1.2) 

n=2 

The Reggeon interaction has interesting properties in the generalized LLA. If one introduces a 
fictitious time in the t— channel, then at the same moment of "time" the interaction occurs only 
between two Reggeons (pair-wise interaction) [7, 8, 9]. It conserves the number of Reggeons 
(elastic scattering) but changes color and two-dimensional transverse momenta of Reggeons. As 
a result, the diagrams contributing to the amplitude An{s,t) have a conserved number n of 
Reggeons in the t— channel and have a form of generalized ladder diagrams [7, 8, 9]. These 
diagrams can be efi'ectively resumed using the Bethe-Salpeter approach [8, 9] and the resulting 
expression for the scattering amplitude An{s,t) can be represented as a sum of contributions of 
the n Reggeon compound states propagating in the channel 

A^{s,t) = tsY^s^-^i^y &,it)PiinAt) . (L3) 
{1} 

Here, En^{q} is the energy of the n Reggeon compound state \Xn,{q}), which satisfies the (2+1)- 
dimensional Schrodinger equation, the so-called Bartels-Kwiecinski-Praszalowicz equation [8, 9] 

'Hn\Xn,{q}) = En,{q}\Xn,{q}) (1-4) 

with {q} being the set of quantum numbers parameterizing all possible solutions. The hamilto- 
nian Tin corresponds to the elastic pair-wise interaction of n Reggeons. It acts on 2-dimensional 

transverse momenta of Reggeons and describes the evolution of the n Reggeon compound state 
\Xn,{q}) i'^ ^he t— channel. The residue factors in (1.3) measure the coupling of the Reggeon 
compound states to the hadronic state and they are defined as 

Ait) - {Mt)\XnM) ■ (1.5) 

For short-distance hadronic states like onium /^^^/^(t) can be calculated in perturbative QCD 
[3, 8, 4]. Expression (1.3) for the scattering amplitude is consistent with the Regge theory 
predictions provided that we interpret the Regge poles as the n Reggeon compound states and 
identify the intercept as their maximal energy 

a{t)-l = max^^jK.M ■ (1-6) 

Combining relations (1.2) and (1.3) together we obtain the Regge asymptotics of the scattering 
amplitude in the generalized LLA as an infinite sum over all possible n Reggeon compound 
states propagating in the t— channel. For n — 2 the first term in the sum (1.2) corresponds to 
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the BFKL Pomeron. The next term, n — 3, is associated with three Reggeon compound states 
which belong to the Odderon family of the Regge poles [10]. Although they have a zero color 
charge, in contrast with the BFKL Pomeron their charge conjugation is negative, C = —1. As 
a result, they cannot couple to the hadronic states with C = 1 like virtual photon in the deep 
inelastic scattering, but for the hadronic states with C = —1 like proton they contribute to the 
total cross sections cr*°* and a^. Moreover, their contribution is responsible for the increasing 
with the energy of the difference [10], Act = — a^. The Regge behaviour of Act is controlled 
by the intercept of the Odderon [11, 12], which despite of a lot of efforts is unknown yet. 

As it follows from (1.2), the contributions of the n Reggeon states to the scattering amplitude 
is suppressed by powers of ctg with respect to that of the BFKL Pomeron. Nevertheless, they 
have to be taken into account in order to preserve unitarity of the yS"— matrix of QCD [7, 8] . To 
derive the higher Reggeon compound states we have to solve the (2+l)-Schrodinger equation 
(1.4) for n interacting Reggeons. For n = 2 it was done in [3], but for n > 3 the problem becomes 
extremely complicated partially due to interaction between color charges of Reggeons. The 
latter interaction can be drastically simplified by taking the multi-color limit [13], ^ oo and 
asN — fixed, and passing by means of Fourier transformation from two-dimensional transverse 
momenta of Reggeons, k± — {ki,k2), to two-dimensional impact parameter space, b± — (61,62)- 
After these transformations, the Reggeon hamiltonian takes the following form [14] 

Hn = ^(i/n + Hn) + 0{N-') . (1.7) 

The hamiltonians Hn and Hn act on holomorphic and antiholomorphic coordinates of the Reggeons 
in the impact parameter space, z — bi + ih2 and z = 61 — ^62, respectively, and describe nearest- 
neighbour interaction between n Reggeons 

n n 

Hn — ^ H{Zk, Zk+l) , Hn = ^ H{Zk, Zk+l) , [Hn, Hn] — , 

fc=l k=l 

where z^ and are holomorphic and antiholomorphic coordinates of the k—th Reggeon and 
periodic boundary conditions Zn+i = Zi and Zn+i = Zi are imposed. Here, H{zi, Z2) +H{zi, z^) is 
the interaction hamiltonian of two Reggeons with coordinates {z\,z\) and (-22,^2) in the impact 
parameter space, the so-called BFKL kernel [5, 14], 

ii{z^, Z2) = -V'(- J12) - ^(1 + J12) + 2^{1) (1.8) 

where il>{x) — ^ Inr(x) and the operator J12 is defined as a solution of the equation 

Ji2(l + J12) = -{zi - Z2fdid2 

with dk — d/dzk- The expression for H{zx, Z2) is similar to (1.8). Thus, in the multi-color limit, 
the 2-dimensional Reggeon hamiltonian 7i„ describing the pair-wise interaction of n Reggeons 
turns out to be equivalent to the sum (1.7) of two 1-dimensional mutually commuting hamilto- 
nians Hn and Hn- This allows us to reduce the original (2+l)-dimensional problem (1.4) to the 
system of two (l+l)-dimensional Schrodinger equations [14] 

Hn\Vn,{q}) = ^n,{q}Wn,{q}) , Hn\<f>n,{q}) = ^n,{q}\<fn,{q}) (1-9) 
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where the wave functions (pn,{q} and <fn,{q} depend only on holomorphic and antiholomorphic 
coordinates of the Reggeons, respectively. Then, in the multi-color limit the spectrum of the n 
Reggeon states can be found as 



Although equations (1.9) are much simpler than the original Schrodinger equation (1.4) it is 
not obvious that they can be solved exactly for an arbitrary number of Reggeons. A significant 
progress has been achieved recently [15, 16] after it was realized that the system (1.9) of (1 + 
1)— dimensional Schrodinger equations has interesting interpretation in terms of one-dimensional 
lattice models [17]. 

Let us consider the one-dimensional lattice with periodic boundary conditions and the number 
of sites, n, equal to the number of Reggeons. We parameterize n sites by holomorphic coordinates 
Zk {k = l,...,n) and introduce the interaction between nearest neighbors on the lattice with 
the holomorphic Reggeon hamiltonian Hn, (1-8). Thus defined "holomorphic" lattice model is 
described by the same Schrodinger equation as in (1.9) and it obeys the following remarkable 
property [15, 16]. It turns out to be equivalent to the celebrated XXX Heisenberg magnet of 
spin s = corresponding to the principal series representation of the noncompact conformal 
SL{2,(C) group. The same is true, of course, for the antiholomorphic Schrodinger equation in 
(1.9). Therefore, the n Reggeon compound states in multi-color QCD share all their properties 
with the eigenstates of the XXX Heisenberg magnet defined on the periodic one-dimensional 
lattice with n sites and their intercept (1.6) is closely related to the ground state energy of the 
magnet.^ This result opens a possibility to apply the powerful methods of exactly solvable models 
[18, 19, 20, 21] to the solution of the Regge problem in QCD. The first step has been undertaken in 
[15, 22] where the generalized Bethe ansatz has been developed for diagonalization of the Reggeon 
hamiltonians. The expressions for the n Reggeon wave functions and their corresponding energies 
were found in terms of the fundamental Q— function which satisfies the Baxter equation for the 
XXX Heisenberg magnet of spin s = — 1. The solution of the Baxter equation was found in the 
simplest case of n = 2 Reggeon state, or BFKL Pomeron. For higher Reggeon states the problem 
becomes much more complicated and it is still open [23]. In the present paper we continue the 
study of the Baxter equation initiated in [15, 22] and develop a method which allows us to find 
its general solutions in the form of asymptotic expansion similar to the well-known quasiclassical 
approximation in quantum mechanics. 

The paper is organized as follows. In Section 2 we summarize the Bethe ansatz solution for n 
Reggeon compound states and introduce the Baxter equation. We interpret the Baxter equation 
as a discrete one-dimensional Schrodinger equation and identify inverse conformal weight of the 
Reggeon states as a small parameter which plays a role of the Planck constant. The quasiclassical 
expansion of the solution of the Baxter equation in powers of this parameter is performed in 
Section 3. Similar to the situation in quantum mechanics it leads to the asymptotic series for the 
energy of Reggeon states whose properties are studied in Section 4. The asymptotic expansions 
for the Reggeon quantum numbers are derived in Section 5. In Section 6 we use the obtained 
results to estimate the energy of higher Reggeon states. Summary and concluding remarks are 
given in Section 7. The analj^ical properties of the energy of the Reggeon states are considered in 

^Indeed, from point of view of lattice models it is more natural to change a sign of the Reggeon hamiltonian, 
Hn, and interpret (—max En) in (1.6) as a ground state energy of the XXX magnet. 
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Appendix A. Relation between Reggeon states and conformal operators is discussed in Appendix 
B. Some useful properties of the asymptotic series are described in Appendix C. 



2. Generalized Bet he Ansatz 

The fact that the system of Schrodinger equations (1.9) is completely integrable implies that 
there exists a family of "hidden" holomorphic and antiholomorphic conserved charges, {q} and 
{q}, which commute with the Reggeon hamiltonian and among themselves. Their explicit form 
can be found using the quantum inverse scattering method as [16, 15] 

Qk — ^ ] ^ ^iii2^i2i3 • • • -^1^11^11^12 ■ ■ ■ ^ifc (2.1) 

n>ii >i2>--->ife>l 

with Zij = Zi — Zj and the expression for qt is similar. The appearance of these operators is 
closely related to the invariance of the Reggeon hamiltonian (1.7) under conformal SL{2,<C) 
transformations [5] 

az + b _ az + b ,^ ^, 

z ; , z — J (2.2) 

cz + a cz + a 

where ac — bd = dc — bd — 1. Indeed, we recognize q2 and q2 as the quadratic Casimir oper- 
ators of the SL{2,<C) group while the remaining conserved charges {qk,qk}, k = 3,...,n can be 
interpreted as higher Casimir operators. The Reggeon compound states belong to the principal 
series representation of the SL{2,€) group and under the conformal transformations (2.2) they 
are transformed as quasiprimary fields with conformal weights {h,h) [24] . 

The n Reggeon states diagonalize the operators {q, q} and the eigenvalues of the conserved 
charges q2, qs, ■ ■ ■, qn play a role of their additional quantum numbers. In particular, the eigen- 
values of the quadratic Casimir operators are related to the conformal weights of the n Reggeon 
state as 

?2 = -h{h - 1) , ^2 = -h(h - 1) , 92 = ?2 ■ 

where the possible values of h and h can be parameterized by integer m and real u 

2 -|- jji _ \ — jYi 

h— — l-iv, h — — l-iv, m — '2. , u — lR. (2.3) 

As to remaining charges, their possible values also become quantized. The explicit form of the 
corresponding quantization conditions is more complicated and will be discussed in Sect. 5. 

To find the explicit form of the eigenstates and eigenvalues of the n Reggeon compound states 
corresponding to a given set of quantum numbers {q, q} we apply the generalized Bethe ansatz 
developed in [15, 22]. The Bethe ansatz for Reggeon states in multi-color QCD is based on the 
solution of the Baxter equation 

A(A)g(A) = (A + i)"Q(A + i) + (A - i)"(5(A - i) . (2.4) 

Here, Q{\) is a real function of the spectral parameter A, A(A) is the eigenvalue of the so-called 
auxiliary transfer matrix for the XXX Heisenberg magnet of spin s = 

A(A) = 2A" + 92A"-' + ... + qn (2.5) 
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and n is the number of Reggeized gluons or, equivalently, the number of sites of the one- 
dimensional spin chain. For a fixed n it is convenient to introduce the function 



g(A) = A"g(A) 



(2.6) 



and rewrite the Baxter equation (2.4) as 

g(A + i) + Q(A -i)- 2Q{X) = f-^^i^ + I + . . . + Q(A) . 



A'- 



(2.7) 



Once we know the function (5(A), the energy of the n— Reggeon compound state can be 
evaluated using the relation 

^^ReeniKqs,. . . ,qn) (2.8) 



where the holomorphic energy is defined as 



dX g A + i) 



(2.9) 



A=0 



The expression for the wave function of the n Reggeon states in terms of the function Q can be 
found in [15, 22] and it will not be discussed in the present paper. 

The Baxter equation (2.7) has the following properties [22]. Wc notice that for g„ = it 
is effectively reduced to a similar equation for the states with n — 1 Reggeized gluons. The 
corresponding solution, g(A), gives rise to the degenerate unnormahzable n Reggeon states with 
the energy 

en{h, qs, qn-i, 0) = £„-i(/i, qs, ?n-i) ■ (2.10) 

These states should be excluded from the spectrum of the n Reggeon hamiltonian and solving 
the Baxter equation for n Reggeized gluons we have to satisfy the condition 

Qn^O- (2.11) 
As a function of the quantum numbers, the holomorphic energy obeys the relations 

Enih, gs, qn) = Enil - h, ^3, g„) = Enih, -qs, (-)"g„) , (2.12) 

which follow from the symmetry of the Baxter equation (2.7) under the replacement h ^ 1 — h 
or A — > —A and qk —>■ {—^qk- This relation means that the spectrum of the Reggeon hamiltonian 
is degenerate with respect to quantum numbers /t, ^3, q^. Then, assuming that the ground 
state of the XXX Heisenberg magnet of spin s = is not degenerate we can identify the quantum 
numbers corresponding to the maximal value of the Reggeon energy as [22] 



h 



1 
2 



(2.13) 



and for the states with only even number of Reggeons, n = 2k, 



Q3 



— ■■■ — Q2k-1 



= 0. 
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For the states with odd number of Reggeons the latter condition is not consistent with (2.11). 

We notice that the conformal weight h enters as a parameter into the Baxter equation (2.4) 
and, in general, one is interesting to find the solutions of (2.4) only for its special values (2.3). 
Moreover, since the equation (2.7) is invariant under the replacement h ^ 1 — h we may restrict 
ourselves to the region 

Re/i>l/2, (2.14) 

or equivalently m > in (2.3). 

Our strategy in solving the Baxter equation will be the following [15, 22]. We will first try 
to solve (2.4) for integer positive values of the conformal weight /i = 2+ and then analytically 
continue the resulting expression for the energy (2.9) to all possible values (2.3) including the 
most interesting one (2.13). It is clear that this procedure is ambiguous since for any integer 
Z one could multiply the energy by a factor exp^inZh). However, according to the Carlson's 
theorem [1], the holomorphic energy is uniquely defined in the region (2.3) by its values at the 
integer positive values of the conformal weight h provided that the function is regular in (2.3) 
and its asymptotics at infinity is 0(e'^^'^^^ with c < tt. As we show in Appendix A, these two 
conditions are indeed satisfied. 

The Baxter equation (2.4) has two linear independent solutions and in order to select only 
one of them we impose the additional condition on the function Q{\) for h — 2+ 

Q(A) A'^-". (2.15) 

For integer positive conformal weight, h > n, the solution, Q(A), of the Baxter equation (2.4) 
under the additional condition (2.15) is given by a polynomial of degree h — n m the spectral 
parameter A. This implies in turn that Q{\) is a polynomial of degree n in A, which can be 
expressed in terms of its roots Ai, A„ as follows [15, 22] 

g(A) = n(A - A,) = exp f^ln(A - A,)) A'' . (2.16) 
j=i \«=i / 

Substituting (2.16) into (2.7) and putting A = A^ we obtain that the roots satisfy the Bethe 
equation for the XXX spin chain of spin s = 

l-n ,'~^^^ k=l.....n. (2.17) 

There is, however, one additional important condition on the roots, which follows from the 
definition (2.6). Comparing (2.16) with (2.6) we find that Ai, \h-n are roots of Q{\) and 
\h-n+i = . . . = A/i = 0. This means that among all the solutions of the Bethe equation (2.17) 
we have to select only those which have the n— time degenerate root A^ = 0. 

The exphcit solution of the Baxter equation (2.4) is known only for n = 2 [15, 22] 

g„=2(A; h) = i^h{l -h),F,(^^ ^' ^ ^ ~ '^ ij (2.18) 

and it can be identified as a continuous symmetric Hahn orthogonal polynomial [25] . Substituting 
the solution (2.18) into (2.9) we obtain the holomorphic energy of n = 2 Reggeon states as [22] 

£2(/i) = -4[V^(/i)-V'(l)] (2.19) 
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where ■0— function was defined in (1.8). We substitute e2{h) into (2.8) and analytically continue 
the result from integer h to all possible complex values (2.3) 



TT 



+ |m| ^ . \ 



This relation coincides with the well-known expression [3] for the energy of the n = 2 Reggeon 
compound state, the BFKL Pomeron. The maximum value of the energy, — In 2, is 

achieved &t h — 1/2 and it is in agreement with (1.1), (1.6) and (2.13). 

Solving the Baxter equation (2.4) for n > 3 one may try to look for the solution as a linear 
combination of the n = 2 solutions [22], Q{\) = J2k Ck{h, q^, .... qn)Qn=2{^\ k). Using the proper- 
ties of Qn=2 as orthogonal polynomials we obtain that the coefficients Ck{h, q^, g„) satisfy the 
multi-term recurrent relations. Most importantly, the quantum numbers qs, ■ ■ ., qn corresponding 
to the polynomial solutions of the Baxter equation turned out to be quantized and their possible 
values as well as the values of the roots were found to be real [22] 

ImAfc = Imgfc = 0. (2.20) 

Although the explicit form of the recurrent relations for {Ck} is known [22] and it is not difficult 
to solve them numerically for lowest values of the conformal weight h and then find the quantized 
charges {qk} and the energy £„, the analytical expression for the solution Q{X) similar to (2.18) 
is still missing. At the same time, the results of numerical solution of the Baxter equation for 
n = 3 and n = 4 presented in [22] indicate that the analytical solution should exist. In this paper 
we will find such kind of solution for an arbitrary n in the special limit, /i ^ 1. We will show 
that it is possible to develop an asymptotic expansion for the solutions of the Baxter equation, 
as well as for the energy of the n Reggeon states, in powers of 1/h. 

2.1. The origin of the quasiclassical approximation 

Let us rescale in (2.7) the spectral parameter as A — > A/i in order to get rid of a large factor 
h{h — 1) in the r.h.s. of the Baxter equation (2.7) as /i — > oo. Then, we introduce new functions 
/(A) and $(A) 

/(A) = exp (/i$(A)) = e"'^^'^'^ Q{hX) , (2.21) 
where according to the definition (2.16) 

*(A)4gln(A-g. (2.22) 

Substituting (2.21) into (2.9) we find the following expression for the holomorphic energy of the 
n Reggeon state in terms of the function $(A) 



(2.23) 



where prime denotes a derivative with respect to the spectral parameter A. 

The main reason for considering the function $(A) is that in the limit /i — > oo we expect the 
scahng Aj = 0{h) and $(A) ~ which allows us to expand $(A) in powers oi 1/h as 

<I>(A) = <l>o(A) + i<l>i(A) + ^$2(A) + ... (2.24) 
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with all functions $i(A) being /t— independent. 

The substitution of the relation (2.21) into (2.7) yields the following equation for /(A) 

/(a + + / (a - - 2/(A) = -V{X)f{X) , (2.25) 
where the notation was introduced for the "potential" 

Taking a naive limit /i — > cxd in (2.25) we find that /(A) satisfies the one-dimcnsional Schrodinger 
equation, — -p^/ = 2m{E — V{X))f, for the wave function of a particle with mass m — 1/2 
and energy E — moving in the singular potential V{X). In this equation 1/h plays a role of the 
Planck constant h and the relation (2.21) has a form of the WKB ansatz, /(A) = exp (^i$(A)). 
This suggests that h ^ oo limit of the Baxter equation (2.25) should be closely related to the 
quasiclassical approximation h ^ in quantum mechanics. 

We would like to notice that there is an intriguing relation between the Bethe Ansatz solution 
of our model and that of the one-dimensional Toda chain [26] . In both cases after the separation of 
variables one has to solve the Baxter equation which has the form of one-dimensional discretized 
Schrodinger equation (2.25) for a particle in an external potential. The only difference between 
models is in the form of the potential and in different boundary conditions which one imposes 
on the solution of the Baxter equation. 

The condition for all terms in the expression for the potential (2.26) to have the same be- 
haviour at large h implies the following scaling of the holomorphic quantum numbers 

|.,3+l,... + _L,» + 0(i), .... |.,„^i,...^J,,,.^0(l). (2.27) 

with % = qf\ Qn = Qn^ and all coefficients q^^^ being /i— independent. This leads to the 
decomposition of the potential V{X) similar to that in (2.24) 

V{X) = Vo{X) + ^V^{X) + ^V^iX) + . . . (2.28) 

where 

Vo{X) = X-^ - q^X-' - . . . - qnX-\ (2.29) 
V,{X) = -A-^-g^^^A-3-...-g«A-", 
V^iX) = -qfX-^-...-q^'X-\ ... 

Finally, to solve (2.25) and find the functions $j(A) wc have to substitute (2.24) into (2.21) and 
(2.25), take into account the decomposition (2.28) and equate the coefficients in front of different 
powers of l//i in the both sides of the discrete Schrodinger equation (2.25). 



3. Leading large— /i approximation 

As the conformal weight h increases the number of roots of the solution of the Baxter equation 
(2.16) and (2.22) also increases. According to (2.20), the roots are real and in the limit /i — > oo 
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it is convenient to introduce their distribution density on the real axis as 



pi<^) = -flY.H^-i^] (3-1) 

with the normahzation condition ^ 

/oo 
da p{a) = 1 . (3.2) 
'OO 

In the large— /i limit p(cr) becomes a smooth positive definite function of a real parameter a. It 
vanishes outside some (not necessary connected) finite interval S on the real axis 

p{a) 7^ , for a e 5 = [(Ti, (Ta] U [0-3, CT4] U . . . U [aafe-i, (T2k] ■ (3.3) 

with the number k depending on the "potential" V{\), or equivalently on the quantum numbers 
^3, Qn- The explicit form of p{a) can be obtained from the Baxter equation (2.25) in the limit 
/i — > 00 as follows. Using (2.22) and (3.1) we get 

$(A) = j^da p(a) ln(A - a) 

and after differentiation of the both sides with respect to A 



where integration is performed along the finite interval (3.3). The last relation implies that ^'{X) 
is an analytical function of the spectral parameter A in the complex A— plane with the cut along 
the finite interval S and singularity at A = 00. Taking discontinuity of $'(A) across the cut we 
can reconstruct the distribution density of roots as 

p((7) = ;^[$'((j + iO) -$'((j-iO)] . (3.5) 
27r 

Similar to (2.24), the distribution density can be expanded in powers of 

P(o-) = Po((7) + -^Pil^) + /^/'2((t) + ... . (3.6) 

Here, the functions Pk{(y) are defined using (3.5) as a discontinuity of the functions ^k{^) and 
they satisfy the normalization conditions, J^dapQ{a) = 1 and J^dapk{<j) = for A; > 1, which 
follow from (3.2). 

To find the solution of the Baxter equation (2.25) in the leading large— /i hmit we use the 
identity 

^ eMh mX ± i/h) - $(A)]) = exp(±i$[,(A) + 0(l/h)) (3.7) 
and keep the leading l//i— terms in the both side of (2.25) to get 

4sin2Q$'o(A)) =yo(A). (3.8) 
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The equation has two different solutions for $o(A) but only one of them satisfies the additional 
condition (3.4) as A — > oo, 



$;(A) = -2i\n 




(3.9) 



Notice that $o(A) can be defined from (3.8) up to a constant 47rZ which does not contribute, 
however, neither to the energy (2.23), nor to the distribution density (3.5). 

Let us consider the properties of $o(A) in the complex A— plane. Taking into account that 
the potential Vo(A) is a real function of A, we find from (3.9) that $o(^) ^'^^ analytical in the 
regions of A, in which the arguments of the square roots change their sign, that is for Vo(A) > 4 
or Vo(A) < 0, and in compact form 

\Vo{X)-2\>2, XeS. (3.10) 

Solving this inequality we find the interval S on the real axis across which the function $o(A) has 
a discontinuity and, as a consequence of (3.5), the distribution density of roots (3.3) is different 
from zero. The same result can be easily deduced from (3.8) as follows. In order to have roots 
on S, the function /(A) should be oscillating on this interval or equivalently $o(A) should have 
an imaginary part for A e 5. If A is outside S, then the function $o(A) is real and it follows from 
(3.8) that < Vo{X) < 4. 

Once we solved (3.10) and obtained the explicit form of S, we substitute (3.9) into (3.5), 
evaluate the discontinuity of '^*o(A) across the interval S and get the distribution density as 



l-lvo{a) + ^-^Vo{a)j , for Vo{a) < (3.11) 




and 



p,(a)^^lnU^Voia)-l + J^Voia)], for Voia) > 4 . (3.12) 



These expressions define the distribution density of roots in the leading large— /i limit for an 
arbitrary number of Reggeons, n >2. 



3.1. Critical values of quantum numbers 

Let us consider the special case n = 2, in which we can compare (2.21) and (3.9) with the exact 
solution (2.18) of the Baxter equation. For n = 2 we substitute Vo(A) = A~^ into (3.10) and find 
the explicit form of 5 as 

po(<7)^0, for S^[-ll]. 

Then, it follows from the definition (3.1) that the zeros of (5n=2(A; h) are located on the interval 
—h/2 < Aj < h/2. Substituting Vo{X) into (3.12) we obtain the distribution density of zeros of 
the n — 2 solutions of the Baxter equation in the leading large— /i limit as 

^•^^"^^^"13^^' -^'^^[-i,!]. (3.13) 
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P(o) 

3.5 - 




1- 

,5- 



-o'.6 -0.4 -0.2 u'o 0.2 0.4 'o. 6 

o 

Figure 1: The histogram of the normahzed roots of the solution of the n — 2 Baxter equation for 
/i = 50 versus the distribution density p{(7) in the leading large— /i limit. 

One can check that this expression satisfies the normalization condition (3.2). To test (3.13) 
we use the explicit expression (2.18) for Qn=2{^',h) at h = 50, find numerical values of all 48 
normalized roots, {\k/h}, and compare on fig. 1 the histogram of their distribution with (3.13). 

Let us generalize our analysis and consider the possible solutions of (3.10) for n = 3 Reggeon 
states. For n = 3 the potential (2.29) is given by Vo(A) = — gsA"^ and, depending on the 
value of the rescaled quantum number ^3, there are two different solutions of (3.10) shown on 
fig. 2. For g| < {q^y the equation Vo{X) — 4 has three real roots A = (Ji,C3,(J4, while for 




(a) (b) 

Figure 2: Possible forms of the potential Vo{X) entering into n — 3 Baxter equation which lead 
to two different solutions of (3.10): (a) ^3 = 0.18 < ^3 and (b) ^3 = 0.25 > q^. The shadow area 
represents the support S. 
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^3 ^ ill) there is only one real root, A = cti. The "critical" value of the quantum number 



27 

can be found as a solution of the system of equations, Vq'(A*) = and Vo(A*) = 4. Notice that 
in both cases, g| < (^3)^ and ^3 > (^3)^, the equation Vo{X) — has only one solution (72 = ^3. 
Then, taking for simplicity ^3 > we identify the finite support S for the distribution density of 
roots at n = 3 as 

5= [(Ji,(72]U[a3,a4], for < ^3 < ^ (3.14) 



with (7i < < (72 = ^3 < (J3 < 1^3 < (T4 and 



(3.15) 



<S = [0-1,(72] , for q3>-^ 

with (7i < < (72 = gs. The generalization of these relations to the case ^3 < is obvious. 
The distribution densities corresponding to (3.14) and (3.15) are shown on fig. 3(a) and (b), 
respectively. 

P(CT) p(CT) 





(a) 



(b) 



Figure 3: Two different forms of the distribution density for the n — 3 Baxter equation: (a) 
?3 < Q3 and (b) g3 > corresponding to the potentials of fig. 2 (a) and (b), respectively. 

It is important to note that, in contrast with the previous case, for n = 3 there exists a 
"critical" value of the quantum number (73 = at which the gap in the distribution of the 

roots vanishes and (3.14) is replaced by (3.15), or equivalently fig. 3(a) is replaced by fig. 3(b). 
One might expect that this value separates two different "phases" of the Baxter equation and 
its solution should have different properties for ^| < ^ and g| > ^. To show that this is really 
the case we recall that we are interesting in polynomial solutions of the Baxter equation (2.16) 
under the additional condition that they should have n— time degenerate root A^ = 0. Till now 
we did not satisfy this condition and "good" solutions could appear in our analysis together with 
"unwanted" ones. In order to get rid of them we have to impose additional constraints on the 
solutions of the Baxter equation in the large h limit. Our assumption is that for n = 3 this can 
be achieved by restricting the possible values of the quantum number ^3 to lie in the region 

1 



<q3< 



1 



27 



(3.16) 
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This condition is equivalent to the statement that for the n — 3 Reggeon compound states the 
distribution density of roots, p(cr), should have the support S given by (3.14) and consisting of 
two connected intervals (see fig. 3(a)). To confirm our assumption (3.16) we present on fig. 4 the 
results (see eq. (6.46) and fig. 8 in [22]) of the numerical solution of the quantization conditions 
for Qs/h^, corresponding to n = 3 and integer conformal weight 3 < /i < 40. We find that all 
numerically found quantized values of qs/h^ are in perfect agreement with (3.16). They lie inside 
the strip defined by (3.16) and for large h their maximal/minimal values asymptotically approach 
the critical values qs/h"^ = Due to the symmetry ^3 —53 of the Baxter equation, (2.12), 

the quantized ^3 are distributed symmetrically with respect to the axis ^3 = 0. Moreover, looking 
on fig. 4 it is difficult do not notice that the possible values of ^3 for different h try to form a 
family of one-parametric curves in the two-dimensional (/i, ^3)— plane which can be labeled by 
integer k — 0,1, ... . An example of such a curve is shown on fig. 4(a). The curves start at 
the points with h = k + 3 and ^3 taking the absolute minimal value, cross the axis ^3 = at 
h = 2k + 3 and for /i — 00 they asymptotically approach the maximal value inside the strip 
(3.16). This suggests that all possible values of the quantum number q^ can be parameterized in 
the following form 

^ = ao{k) + ^a,ik) + ^a2{k) + 0{h-') , A; = 0, 1, ... (3.17) 

where aj{k) are some unknown coefficients depending on integer k. Comparing this ansatz with 
(2.27) we find that ao(A;) = ^3 and all remaining coefficients ai{k), a2{k), ... correspond to higher 
l//i— terms in the expansion (2.27). To test (3.17) we have to evaluate the nonleading 1/h 
corrections to the quantum number qs/h^ in (2.27). This will be done in Sect. 5. 




Figure 4: The distribution of quantized q^ for different values of the conformal weight 3 < /i < 40: 
(a) The results of the numerical solution of the n — 3 Baxter equation. The dashed fine represents 
the upper and lower limits for q^ in (3.16). The solid line indicates one of the curves to which 
quantized ^3 belong, (b) One-parametric families of curves defined in (5.15) and (5.16). 

Let us generalize the condition (3.16) for the Baxter equation at n > 3. For an arbitrary n 
the solutions of (3.10) have the following properties. Since Vo(A) ~ 1/A^ — 0+ as A — cxo the 
interval S does not contain infinity and has a finite size. Then, the potential is singular at the 



15 



origin, Vo(A) ~ —Qn/^"' — > oo as A — > 0, and, according to (3.10), the point cr = is always 
inside S. In general, depending on the value of the quantum numbers gs, the interval 

S may consist of up to n — 1 connected intervals. We require that the number of connected 
intervals inside S should take a maximal possible value, n — 1, for the polynomial solutions of the 
Baxter equation.^ The necessary condition for this is that the equations Vo(A) = and Vo(A) = 4 
should have n — 2 and n real roots, respectively. The solution of these constraints leads to the 
quantization conditions on q^, . . ., qn similar to (3.16). To satisfy them the quantized values of 
the conserved charges should lie inside the domain in the n — 2 dimensional space of all possible 
^3, . . ., g„. The boundary of the domain defines the "critical" n — 3 dimensional hypersurface 

Sn(93,---,?:) = (3.18) 

which separates different phases of the Baxter equation. For n = 3 we have ^3(^3) = (q^)'^ — 27 
and the quantized q^ belong to the interval (3.16). Let us find the explicit form of the quantization 
conditions for n = 4. 

For n — A the potential is given by Vo{X) — A~^ — qs — ^4 A"^ and in order for the equation 
Vo(A) = to have two real roots we have to require that ^4 > — Under this condition the 
function Vo(A) has two extremums - two local maximums for ^4 > and one maximum and one 
minimum for ^4 < 0. Their positions, A = A 

max ov A — Amin) Can bc fouud as solutions of the 
equation ^ = 0. Then, one can easily check that the second condition on Vo(A), namely that 
Vo{X) — 4 has four real solutions, becomes equivalent to the requirement Vb(Ainax) > 4. After 
some algebra we find that the second constraint on ^3 and q^ can be represented together with 
the first one in the following form 

These inequalities define the two-dimensional domain on the (^3, ^4)— plane shown on fig. 5. 
Using the results of numerical solutions of the n = 4 Baxter equation [22] we plot on fig. 5 the 
numerical values of quantized {q'i/h? ,qi/h*) corresponding to the conformal weight h — (see 
fig. 10 in [22]). We observe a complete agreement of the numerical results with (3.19). Moreover, 
the quantized values of ^3 and ^4 arc distributed on fig. 5 in a regular way, but similar to the 
previous case n = 3, in order to describe their "fine structure" we need to know nonleading 1/h 
corrections to both quantum numbers. 

The one-dimensional boundary of the domain (3.19) defines the "critical" values of the quan- 
tum numbers. The range of quantized charges ^3 and ^4 can be found from (3.19) as 

1.1 1.1 

<03<— , <0'4< . 

4 - - 4 ' 64-^^-16 

We notice that taking ^4 = in (3.19) we obtain the relation g| < 1/27, which wc identify as the 
condition (3.16) on ^3 for n = 3. In the same manner, the point 53 = 54 = on fig. 5 corresponds 
to the possible range of the quantum numbers for the n — 2 Baxter equation. Continuing this 

^This condition has its analog in the Toda chain [27, 28]. To separate the variables and reduce the Bethe 
Ansatz to the solution of the Baxter equation, one performs the canonical transformation and introduces the set 
of n — 1 generalized coordinates X\, Xn-i [29]. Then, the n — 1 closed intervals inside S correspond to the 
regions of the classical motion of x\, Xn-i- 
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-0.2 -0.1 0.1 0.2 

Figure 5: The distribution of the qiiantized ga and corresponding to the polynomial solutions 
of the n = 4 Baxter equation. The quantum numbers lie inside the generalized triangle defined 
in (3.19). 



analogy, we may conclude that the n — 2 dimensional domain of the quantized charges gs, . . ., g„ 
corresponding to the n Reggeon states lies on an intersection of the hyperplane qn+i = with the 
n — 1 dimensional domain of gs, . . ., qn+i corresponding to the (n + 1) Reggeon states. The 
same property can be expressed in terms of (3.18) as Tlniql, = ^n-i-i(?3; • • • , (In+i = 0) • 

3.2. Divergences of quasiclassical approximation 

Let us use the quasiclassical expansion (2.24) of the solution of the Baxter equation to evaluate 
the holomorphic energy (2.23). In order to apply (2.23) we need to know the behaviour of $'(A) 
at A = ±i/h. In the limit /i — > oo, the points A = ±i/h move toward the origin from different 
sides of the imaginary axis and the energy (2.23) is different from zero because the point A = 
belongs to the interval S. Using (3.9) and (2.29) we find the asymptotic behaviour of the function 
$o(A) for A — > at the both sides of the cut as 

%iX±iO)'^'Tiln^ + ... (3.20) 

and applying (3.5) we conclude that the density of roots is logarithmically divergent at the origin 
(see figs. 1 and 3) 

pia)''^' - In ^ + 0{a). (3.21) 

Substituting (3.20) into (2.23) we obtain the following expression for the holomorphic energy 

= -21n(g„/i") ^ -21ng„ , as h ^ oo . (3.22) 

Thus, the holomorphic energy of the n Reggeon compound states scales in the limit h ^ oo 
as a logarithm of the "higher" quantum number g„.^ We stress that (3.22) was found in the 

•^A somewhat similar result was obtained in [14]. 
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leading large— /i limit and in order to justify this expression we have to estimate the level at 
which next-to-leading corrections to the energy appear. Naively one might expect that 
nonleading corrections to e„ will arise at the 0(1/ h) level. However, as we will show in a 
moment, this is not true due to singular behaviour of the potential Vo(A) at A = and, as a 
consequence, the quasiclassical expansion of the energy does not work properly. Indeed, beyond 
the leading large— /i approximation the energy (2.23) gets contribution from $i, $2, • • •, which 
in the quasiclassical approximation are proportional to the derivatives of the function $o(A) (see 
eq. (4.1) below). In particular, ~ '^dx^+i ~ -^"^ as A — and the contribution of to the 
holomorphic energy can be estimated using (2.23) as e„ ~ \~^h~^ = 0{h^) for A = ±i/h and it 
is of the same order as the leading order contribution (3.22). Thus, in order to get an expansion 
of the holomorphic energy in powers of l//i it is not enough to restrict ourselves by the leading 
order result (3.22). 



4. Beyond the leading order 

To find the solution of the Baxter equation (2.25) beyond the leading 1/h order we have to keep 
nonleading terms $1, $2, ... in (2.24) and expand further the r.h.s. of (3.7) in powers of 1/h. 
Substituting the expansion for /(A± ^)//(A) into the Baxter equation (2.25) and comparing the 
coefficients in front of powers of 1/h in the both sides of the equation we obtain the following 
system of equations 

<l>i(A) = g*o+^*!.'*'i + ™t*i(5(»S)'-^(»'.)'-^»'0+ill5r. - (4.1) 

where prime denotes a derivative with respect to the spectral parameter A, the function $o(A) 
was defined in (3.8) and (3.9) and the potentials Vi, V2, ... were introduced in (2.28) and (2.29). 
Integrating these relations with the boundary conditions $jfc(oo) = one can get the functions 
$i(A), $2(A), ... and then reconstruct the solution of the Baxter equation using (2.21) and (2.24). 

Since we would like to use the solution of (4.1) to find the energy of the n Reggeon states (2.23) 
it is of most interest for us to consider the behaviour of $a;(A) at small A = 0{h~^). However, 
as it follows from (4.1) and (3.20), the expansion for the function $(A) in powers of l//i is not 
well defined for small A due to singular behaviour of the potential at the origin, ^(A) ~ 1/A" 
as A — > 0. Indeed, as was anticipated at the end of previous section, the functions ^^(A) are 
proportional to the higher derivatives of $o(A)- For small A one can use (3.20) to find from (4.1) 
the following general form of their asymptotic behaviour as A — > 

$o(A) = a_i ln(a) + oo + oiiA + a2(iA)^ + O(A^) 

$;(A) = bo/iiX)+h + b2iX + 0{X^) (4.2) 

*2(A) = co/{tXy + c,/{tX)+C2 + 0{X), 

wliere the coefficients depend on whether A approaches the origin from above or below the cut 

5. To calculate the holomorphic energy (2.23) we substitute A = ±i/h into (4.2) and evaluate 
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$'(A) using (2.24). In particular, $'(— has the following expansion 

^\-i/h) = -a_i In /i+ (ao + 60 + Co + ■ ■ ■) + («i + ^1 + ci + • ■ •)/^+ («2 + &2 + C2 + . . ■) /h^ + 0{l/h^) 

(4.3) 

where dots denote the contributions of the functions with A; > 3. We conclude that different 
terms in the large— /i expansion (2.24) give the contributions to $'(±i//i) of the same order in 
1/h. As a consequence, in order to calculate the holomorphic energy e„ we have to find an 
effective way to resum the contributions of ^^(A) at small A = 0{l/h) to all orders k. 



4.1. Leading order resummation 

Let us try to simplify the system of equations (4.1) in the limit of small A. As was shown in 
Sect. 3, the function $o(-^) ^ '^^^ along the interval S of the real axis and for A ^ we have 
to consider separately two cases when A approaches the origin from above and below the cut, 
ImA > and ImA < 0, respectively. Using (3.20) we obtain the following identities 

rfl.'n^l*-»/^(l + '^W' for ImA >0 
■^Pf'*"*^" = |£(1 + 0(A)), tor ImA <0 

Without lack of generality we may concentrate on the first case, Im A > 0, and apply the relation 
(4.4) in the form exp(i$Q) ^ exp(— z^q) to simplify the trigonometric functions entering into the 
system (4.1) as 

sm{%) = lexp(^$'o)(l + 0{X'^)) , cot{^',) = i + ^(A^") . (4.5) 

Then, in the hmit A — > and ImA > the system of equations (4.1) is replaced by 

$;(A) = + ^yl(A)e-**o + o(A2«-l), (4.6) 

where $q can be defined from (3.8) as exp(i$Q) = 2 — Vo(A) + (9(A^"). One easily checks 
that performed approximation introduces ambiguities into the definition of the functions $fe(A) 
at the level of C(A^"'~'') corrections, ~ A^""^ ~ A^"~^, .... Therefore substituting 
these relations into (2.23) and (2.24) and putting A = ±i//i one can find the expression for the 
holomorphic energy with the following accuracy 

This means that using the approximation (4.5) and (4.6) one can define the holomorphic energy 
up to e)(/i-2") terms. 

The reason why we are interesting in the approximation (4.5) is that instead of solving the 
infinite set of equations (4.6) for one can find a closed expression for the function $'(A). Our 
observation is that for small A the Baxter equation (2.25) can be effectively replaced by a new 
equation 

"^^Mx'/^^ = 2 - V{X) , for A^OandImA>0 (4.7) 
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and 



/(A) 



2-V{\) 



for A ^ and Im A < , 



(4. 



which reproduces the relations (4.6) and which can be solved exactly. A simplest way to show 
this is to use (4.4) and (3.7) and to take into account that for small A in the Baxter equation 
(2.25) we have f{X-i/h) < f{X + i/h) for ImA > and f{X-i/h) > f{X + i/h) for ImA < 0.^ 
Let us consider the relation (4.7) and use the ansatz (2.21) to rewrite it as 



exp 



/i(et^- l)*(A)] ^2-V{X). 



where d = S^. Its solution has a form 

aX 



i^'{X) = 



exp(f a) - 1 



ln(2-y(A)), for ImA >0. 



(4.9) 



As a check, we substitute $ = J2k>oh~'^^k and V = J2k>oh~''Vk into the both sides of (4.9), 
compare the coefficients in front of powers oi 1/h and reproduce (4.6). An analog of (4.9) for 
ImA < can be obtained from (4.8) as 



i$'(A) = 



exp(-f 9) - 1 



ln(2-1/(A)), for ImA <0. 



(4.10) 



Substituting (4.9) and (4.10) into (2.23) and putting A = ±.i/h we derive after some algebra the 
following expression for the holomorphic energy of the n Reggeon state 



d 



exp ((9) 



In 



2-V\ 



h 



V 



' h 



+ C>(/i-'") 



(4.11) 



A=l 



and taking into account that the potential (2.28) is a real function of A 

d 



-2 Re 



exp (9) 



In 2- V 



tX 
h 



(4.12) 



A=l 



Here, the last term was added to indicate explicitly that the expression for the holomorphic 
energy was found within the approximation (4.5) and it is valid only up to 0(/i~^") terms. This 
implies in particular, that expanding the operator 



d 



exp (9) 



k=0 



k\ 



with Bk being Bernoulli numbers, we may keep in (4.12) only 2n — 1 first terms. 

The expression for the energy (4.12) is real. Another interesting property of (4.12) is that 
although it has a rather formal form it can be explicitly evaluated for an arbitrary potential 
V{X). The calculation is based on the relation 



d 1 

_ In(aA) = In(aA) - — 

ex])(J) — 1 2A 



1 



+ 



1 



12A2 120A4 



+ ...==^(A)+lna 



(4.13) 



''Since in both cases the points A ± i/h should Ue on the same side from the real axis, A cannot approach the 
origin too close, |ImA| > 1/h. 
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which is vahd for an arbitrary a and which follows from the definition of the ■0— function given 
in (1.8) and from identity (In A)' = (e^ — i)ip{X) — ip{X + 1) — — 1/A. Let us consider the 
equation V{X) — 2 with the potential defined in (2.26). It has n roots which we denote as Si, 

Substituting this relation into (4.11) and applying the identity (4.13) we find 

n 

en^-2\n2-Y,['^{'^ + ih5i)+ip{l-ihSi)-2ilj{l)] + 0{h-^''). (4.15) 

i=l 

Notice that the roots 5j depend on h and in the limit /t — > oo they scale as 6 ^ h^. For example, 
for n — 2 one can easily obtain from (4.14) their explicit form as 5i l = 

It is interesting to note that the parameters Si have a simple interpretation in terms of the 
auxiliary transfer matrix A(A) defined in (2.5). Comparing the definitions (2.5) and (2.26) we 
conclude that hSi, hSn entering into (4.15) can be determined as zeros of the auxiliary transfer 
matrix 

A(Mfc) = (/i5fe)"(2- 1^(4)) = 0, A; - l,...,n. 

Using the well-known asymptotic expansion (4.13) of the ■0— function we get from (4.15) the 
expansion of the energy in powers oil/h 

er.--Wn + - 4 E ^^^^^^^^p^'' + 0{h-'n (4.16) 

where ■0(1) = — 7e — —0.577216 is the Euler constant and C{2k) is the Riemann zeta function. 
The comparison of (4.16) and (3.22) shows that, in accordance with our expectations (4.3), 
nonleading l//i— corrections to the solution of the Baxter equation generate 0(/i°) contribution 
to the energy, 2n-0(l), as well as an asymptotic series in 1/h. We should stress that applying 
(4.16) we have to take also into account the /i— dependence of the roots S^ and expand J2^=i S^'^^ 
in powers of 1/h. Although (4.16) looks like an even function ol 1/h the latter expansion will 
give rise to odd powers of 1/h. 

Let us apply (4.16) for n = 2. Using the explicit form of the roots, S^^ — together with 
q'i — —h{h — 1) we obtain from (4.16) 

e,{h) = -2 \n{h{h - 1)) - 47, - ^^^^ _ + 0{h-') = -4 In - 47, + ^ + ^ + 0{h-') . (4.17) 

Comparing this expression with the exact result (2.19) we verify that up to 0{h~'^) corrections 
both expressions do coincide. 

Expression (4.16) defines only 2n — 1 first terms of the infinite asymptotic series for and 
in order to apply (4.16) it is very important to understand how accurately these few terms 
extrapolate the exact result for for an arbitrary value of the conformal weight, 1/2 < h < 00. 
The answer to this question is based on the properties of the asymptotic series [30, 31, 32] and 
it will be discussed in Sect. 6. To anticipate our conclusions, we substitute the most interesting 
value of the conformal weight, (2.13), into (4.17) and (2.19) and obtain the following numerical 
values 

£2''''(l/2) = 5.797059, £2(1/2) = 8 In 2 = 5.545177 . (4.18) 

Remarkably enough, the approximate result turns out to be only 4% bigger the exact expression, 
which determines the intercept of the BFKL Pomeron (1.1). 
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4.2. Improved approximation 



The approximation performed in the previous section has a strong hmitation. It does not allow 
us to evaluate the holomorphic energy with the accuracy better that 0{h~'^"'). Since we do not 
know in advance how many terms in the expansion of the energy we will need in order to 
approach the physically interesting value £n(/i = 1/2) it is desirable to have an improved scheme, 
which would allow us to predict higher 1/h corrections to the energy. 

There is a simple way how one can improve the approximation. Let us start with the Baxter 
equation (2.25) and rewrite it in the following form 



/(A + t/h) 
/(A) 



2 - V^(A) - 



/(A) 



7(A-^A). 



(4.19) 



As we have seen in Sect. 4.1, in the limit of small A and ImA > the ratio /(A + i/h)/f{\) 
behaves as ~ —V{X) ~ qn^^^^- This fact allows us to neglect the last term in the r.h.s. of the 
Baxter equation (4.19) and reproduce (4.7). The next step will be to consider the same term as 
a small perturbation and iterate the Baxter equation (4.19) as 



f{\ + i/h) 



2 - V{\) - 



1 



/(A) 2-V(X-^|h)-{ji^) 
Repeating this procedure A;— times we obtain 



-1 ■ 



f{X + i/h) 



/(A) 

2-V{\-i) 



= 2 - V{X) 

2-V{X-f) 



(4.20) 



^ '^V'^ T) [f(X-i{k+l)/h)) 



1 -1 



1 -1 



Here, the ratio (/(A — ik/h)/ f{X — i{k + l)/h))''^ is much smaller than V{X — ik/h) at small 
A ~ 1/h provided that ImA > and the arguments of / still lie in the upper half plane. 
Otherwise, for ImA < 0, the same ratio becomes large and we are not allowed to treat it as a 
small parameter. This means, that for fixed number of iterations, we have to replace Im A > by 
a stronger condition ImA > {k-\-l)/h. In particular, we cannot continue the iterations infinitely 
and replace the r.h.s. of (4.20) by continued fraction since this will push the allowed region of A 
away from the origin. 

Let us consider (4.20) in the case when A = 0{l/h) and ImA > {k + l)/h. Then, using the 
relations V{X-ik/h) = C(/i") and /(A - ik/h)/f{X - i{k + l)/h) ~ -V {X-i{k+l)/h) = C(/i") 
we can find that neglecting (/(A — ik/h)/f{X — i{k + l)//i))"^ in (4.20) we modify the r.h.s. of 
the Baxter equation (4.20) at the level of C(/i'~^^*^'*"^^") terms. Performing this transformation 
and introducing notation for the "iterated" potential 



Vfe(A) = 2 - V^(A) 



2-V{X-i)- 2-V{X-f) 



2-V{X 



ik 
h 



1 -1 



(4.21) 



we can rewrite the Baxter equation (4.20) in the following form 

f{X + z/h) 



/(A) 



Vfc(A) + 0{h-^^^+'^^) , 
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where ImA > (A; + l)/h and A = 0{l/h). For k — Q this equation is equivalent to (4.7) and 
similar to (4.9) its solution can be represented for an arbitrary k as 

i$'(A) = _ ^ In (V.(A) + 0{h-^-'^^>)) = _ ^ In V.(A) + 0(/^-(^'=+^)-) (4.22) 

where we took into account that Vfe(A) = 0{h^) for A = (9(l//i). We recall that (4.22) was derived 
under the additional conditions IniA > {k + l)/h and A = 0{l/h). To find the expression for 
i$'(A) in the lower half-plane, Im A < 0, we have to iterate the Baxter equation (4.19) considering 
/(A + ^)//(A) as a small parameter and follow the same steps which led to (4.10). The final 
expression is similar to (4.10) with 2 — V^(A) replaced by complex conjugated iterated potential 
Vfe(A) and it is defined for 1mA < -(A; + l)/h and A = 0{l/h). 

Although (4.22) was found for 1mA > {k + l)/h we may analytically continue the result to 
\ = i/h and substitute it into (2.23) to get the following expression for the holomorphic energy 



d 



exp (9) 

and using the fact that V{\) is a real function of A 



-2(fe+l)n 



) 



A=l 



-2 Re 



d 



exp (5) 



+ 0{h 



-2{k+l)n 



)■ 



(4.23) 



A=l 



Here, /c is a positive integer which enters into the definition (4.21) of the potential Vk- We 
conclude from (4.23) that the accuracy of the approximation is controlled by the number of 
iterations k. Increasing this number in the definition (4.21) of the potential and using (4.23) 
we can obtain the expansion of the holomorphic energy in powers of with an arbitrary 
accuracy. For A; = the expressions for the energy (4.23) and (4.12) coincide. Each new iteration, 
A; — >■ /c + 1, adds 2n additional terms to the expansion of e„. 

One can apply (4.23) in a two different ways. The simplest one is to expand lnVyk(iA//i) and 
d/ (exp((9) — 1) in powers oil/h and 5, respectively, and keep only first 2(A; + l)n — 1 terms. The 
second way is based on the identity (4.13). Similar to (4.14) one has to decompose Vjt(A) into 
simple factors and obtain the expression for as a sum of ■0— functions. 



5. Fine structure of quantum numbers 

The expressions for the energy, (4.23) and (4.21), depend on the potential V^(A). Applying (4.23) 
one finds that the higher terms in the asymptotic expansion of the holomorphic energy become 
sensitive to the nonleading 1/h corrections to the quantum numbers qk/hf' entering into the 
definition (2.28) and (2.29) of V{\). The results of the numerical solution of the Baxter equation 
indicate that the nonleading corrections to the quantum numbers are organized in such a way 
that for different values of the conformal weight h the quantized values of ^3, g„ belong to the 
family of one-parametric curves for n = 3 (see fig. 4), two-parametric surfaces for n = 4 and, in 
general, to n — 2 parametric (n — 2)— dimensional hypersurfaces for an arbitrary n. There is a 
simple way how one can understand this fine structure within the quasiclassical approximation. 

Let us recall that in the large h limit, the distribution density of roots of the polynomial 
solutions of the Baxter equation has the support S consisting of n — 1 connected intervals. The 
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total number of roots is equal to the conformal weight h and, as a consequence, the distribution 
density (3.1) satisfies the normalization condition (3.2). All h roots take real values and they are 
distributed on the interval S = [cti, (72] U ... U [<72n-3, o'2n-2]- Let us denote by A^^ the number of 
roots belonging to the interval [(J2k-i,<^2k] with k = 1, n — 1. We remember that one of the 
intervals, say [o'2n-3,c!'2n-2], necessary contains the origin a = and the corresponding root of 
the polynomial solution should be n— time degenerate. This leads to the following conditions on 

n-1 

J2Nk = h, 0<Ni,N2,...,Nn-2<h-n, n < N^-i < h . 

k=l 

We recognize that for a given n and h among all integer numbers A^i, A^n-i there are only n — 2 
independent. It is natural to expect that these are the sets of the integer numbers {A^fc}, which 
parameterize (n — 2)— dimensional hypersurfaces describing the distribution of the quantized ^3, 
Qn for different values of h. 

To show this we use the definition (3.1) of the distribution density and express the numbers 
Nk as 

/ da p{a) = — ^ . 

In the large h limit, we substitute (3.6) into this relation and obtain expansion of its l.h.s. in 
powers of 1/h. At the same time, depending on the value of A^^, the r.h.s. contains only ~ /i° 
and ~ h"^ terms. Then, from the comparison of the both sides we get 

/ da [po{a) + -pi{a)] = — , / da pj{a) = , (5.1) 

where k — 1, ...,n — 1 numerates the intervals inside S and j — 2, 3, ... refers to the level of 
nonleading I//1— terms in the expansion of the distribution density (3.6). We notice that for 
small values of A^^, such that A^^ = 0{h^), the first equation in (5.1) can be further split into two 
independent relations 

/ dapo{a)^0, / dap,{a)^Nk, iorNk^O{h^). (5.2) 

Finally, using the definition (3.5) of the distributions pj as discontinuity of the function across 

the cut, we rewrite the integral of pj over the cut [a2k-i-i C2fc] as a contour integral of around 
the cut 

£;^rf.ftW=4^*;(A) (5.3) 

where the contour Ck encircles the interval [o"2fc-i, o'2fe] in the complex A— plane in the anti- 
clockwise direction. Using the relations (4.1) and replacing the potentials Vj(A) by their defi- 
nitions (2.28) and (2.29), we can express the r.h.s. of (5.3) in terms of quantized qn\ 
Combining (5.1), (5.2) and (5.3) together we obtain an infinite set of equations, in which we treat 
the numbers A^i, A„_i as parameters and ^3', qn^ as unknown variables. For a given set 
of integers A^i, A^n^i their solution will give us the expressions for ^3', qll^ which we will 
substitute into (2.27) to find the asymptotic expansion of the quantum numbers 53, qn- This 
result can be expressed in the following form 

n-1 

/i=^A,-, qk^qk{h;Ni,N2,...,N^-2), Ni, N^_2 > , N„_i > n (5.4) 
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where k — 3, n. If we relax the condition for Nk to be integer, then these parametric 
relations define the n — 2 dimensional domain in the space of quantum numbers /i, 53, g^. 
The interpretation of (5.4) within the framework of two-dimensional conformal field theories is 
proposed in Appendix B. 



5.1. Nonleading corrections at n = 3 

Let us find the explicit form of the function qk = qk{h; Ni, Nn-2) for the n = 3 Rcggeon 
states. For n = 3 the distribution density of roots on S = [o"i,o"2] U [o"3,(T4] is described by two 
integer numbers < Ni < h — 3 and N2 = h — Ni > 3. According to our notations, N2 > 3 
counts the number of roots on the interval [cxa, a^] including the 3-time degenerate root = 0. 
For integer conformal weight h — Ni + N2 one can form only h — 2 possible sets of integers 
{Ni, N2) = {Ni, h — Ni) with iVi = 0, ...,h — 3. This means that in agreement with the numerical 
results shown on fig. 4, for fixed conformal weight h >3 there are only h — 2 possible quantized 
values of q^. 

Let us analyze the system of equations (5.1) in the special limit A^i = 0{h^), in which (5.2) 
holds. The first equation in (5.2) involves the leading distribution density po(c), which according 
to (3.12) and (3.11) is a positive definite smooth function on S (see fig. 3(a)). Therefore, in order 
to satisfy J^^ da po{a) = we have to require that the interval of the integration is shrinking to 
a point, (Ti = (T2. As we have shown in Sect. 3.1, this is the same condition which defines the 
critical value of the quantum numbers. We conclude that the solution of the first equation in 
(5.2) is that the leading quantum numbers ga, should take their critical values qk — qIj 

equivalently belong to the critical hypersurface E„, (3.18). For n — 3 this amounts to say that 

q, = qP^±-^. (5.5) 



Let us consider the second equation in (5.2). It contains the integral of the first non-leading 
correction to the distribution density, pi(cr), over the vanishing interval. In order for the integral 
to be different from zero, the function pi{a) should have singularities on [(71,(72]. To show this 
we introduce a small parameter e, 

a2-ai = e, Vo(ai) = ^0(^2) = 4 (5.6) 

and approach the limit a2 — Ci ^ as e — 0. Here, the second relation follows from (3.10) 
(see also fig. 2(a)) and the potential Vq for n = 3 was defined in (2.29). Choosing for simplicity 
qs > 0, we use (5.6) to obtain the following expansions at small s 

""'''^Tll^l' ^'^72f"V'' Vo{a)=A + lUx{l-x)e', (5.7) 



where a — ai + ex and < x < 1. Substituting these expressions into (3.8) and (4.1) we find the 
asymptotic behaviour of $o(A) and $i(A) at the vicinity of the infinitesimal interval [(7i, (72] as 

$'o(A) = 12e^ z{z - 1) + 0{e') , $;(A) = (; + 7^) + /'/^^ , + C?(6°) (5.8) 
V AeKz z-lj 24eJz{z-l) 
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where \ — ai + ze, z has an infinitesimal imaginary part and < Rez < 1. Taking discontinuity 
(3.5) across the cut < 2; < 1 we obtain the first nonleading correction to the distribution density 



pi{ai + xe) = - 

e 



-i(.(.)..,.-x)).^-^- 



+ C»(£°) . (5.9) 

In accordance with our expectations it is divergent for £ — > 0. Applying the identity 

I da Pj{cr) — £ dx Pj{(7i + xe) — si — — $^((7i + ez) (5.10) 

with Ci encircling the interval [0, 1] in complex 2;— plane, we may use (5.9) to evaluate the integral 
in (5.10). However, a more effective way to obtain the same result is to consider the contour 
integral of ^[{ai+ze) in (5.10) and reahze, using (5.8), that the contour Ci can be deformed away 
from the interval [0, 1]. Then, the contour integral in (5.10) is given by a residue of ^[{ai + ze) 
at 2; = 00 (with ze = fixed!) and its substitution into the second equation (5.2) yields 

Finally, using the definition (2.29) of the potential Vi(A) we obtain the first nonleading correction 
to the quantum number as 

- -— (5.11) 

where A^i = 0{h^) is a positive integer parameterizing all possible solutions of the quantization 
conditions (5.1) and (5.2). 

Let us repeat similar calculation and find the next correction, q^^ . We start with the equation 
(4.1) for $2 use the asymptotics (5.8) to estimate the leading singularity of ^2{<Ji + ze) in 
the limit e — * as ~ e~^. This seems to imply that the integral over P2(c) in (5.10) should 
be divergent at £ — > 0. However, despite of the fact that all these singular terms potentially 
contribute to p2(o'), the integral J^^ dap2{cr) gets a nonzero contribution only from 0{e~^) term 
and is finite. To understand this property we consider the contour integral of $2((7i + ez) 
in (5.10) and take into account, using (5.8), that to any finite order of the expansion of the 
functions $^(o"i + ez) {j = 0, 1, ..) in powers of e their singularities are located on the interval 
[0, 1] and at 2; = 00 in the complex 2;— plane. As a result, the contour integral in (5.10) is given 
by the residue of the function $j (cri + ez) at infinity ^ = 00 but with ez = fixed. One can find 
from (5.8) that in the hmit z — 00 and ez = fixed the functions $j((Ti + ez) have the following 
asymptotic behaviour ~ J2k{£z)~'^ fk with some coefficients fk- All (52;)"'^— terms with k ^ 1 
in the expansion of $2(0"! + ez) have a zero residue at infinity and the contour integral (5.10) 
gets a nonzero contribution only from 0{{ez)~^) term. Carefully expanding (4.1) in powers of e 
we identify after some algebra the proper terms in $2(0-1 + ez) as 

$'2(<7, + £.) {ez)-'(2{q^^r - - I + + ... (5.12) 

where dots denote terms with another (smaller or larger) powers of 1/(^2;). Here, the potential 
V2(A) = —\~^q^^ was defined in (2.29) and ^3^' is the next nonleading correction to the quantum 
number q^. The quantization condition (5.1) for p2 has a form J^^ da p2{cr) = §dz^2 = and 
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it requires that the residue of (5.12) at z — oo should be equal to zero. As a result, the second 
correction to the quantum number can be found from (5.12) as 

and after substitution of (5.11) the explicit dependence of ^3^' on integer number Ni is 



,(2) ^ , 

~ 27 



= — (qn^ + 15Ni + 4) . (5.13) 



Continuing the same procedure it is straightforward to obtain the next nonleading corrections 
to ^3 from the quantization conditions (5.1). Here, we present the results of our computations 
using Maple V Symbolic Computation System 

/3 

qf = -^(2N^ + 12N^ + 32Ni + 8) 
2187 

^ (114 + 330 - 5478 - 10050 - 17012 A^i - 7360) (5.14) 



(30 + 60 - 942 - 972 iVi - 688) 



6561 

9f = fl344 + 5058 - 96210 Nf; - 254382 - 709062 N?; 

^ 59049 V ^ ^ ^ ^ ^ 



-668508A^i- 261472) 

(u 

531441 



/3 

g« = y f 16560 A^7 + 76104 - 1708338 iVf - 5887110 - 23283058 N'f 



-34667562 Nf - 30609348 A^i - 10272032] 

Finally, combining (5.5), (5.11), (5.13) and (5.14) together we obtain the first eight terms of the 
asymptotic expansion (2.27) of the quantized ^3 in powers of h. We notice that (5.11), (5.13) 
and (5.14) were found for positive ^3"'. To get the corresponding expressions for negative ^3°^ 
we may use the symmetry of the Baxter equation (2.12) under ^3 — > — ^3 to change a sign of all 
nonleading coefficients ^3^ in (5.11), (5.13) and (5.14). 



5.2. Comparison with numerical calculations 

Substituting (5.14) into (2.27) we find that the resulting expression for quantized q^ = q3{h;Ni) 
has a parametric form (5.4) with h — Ni + N2, Ni > and N2 > 3. Choosing different values of 
the integers A^i and A^2 we get the spectrum of quantized 53 and h which should be compared with 
the numerical results shown on fig. 4(a). To perform the comparison we remove the condition 
for A^i and N2 = h — Ni to be integer and consider the dependence of ^3 on A^i and h in two 
different cases. In the first case, to which we will refer as to A, we take Ni to be zero or positive 
integer and A^2 > 3 to be positive real. This gives us two families of curves 

A±: qs = ±q3ih;N,), iVi = Z+ , /i > 3 (5.15) 

where ± refers to the symmetry ^3 —q^ and h is a continuous. In the second case, B, we 
choose N2 to be positive integer and Ni = h — N2 to be positive real, 

B± : q3 = ±q3{h;h-N2), = 3 + Z+ , /i > iVs ■ (5.16) 
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As a result, we obtain the families A± and B± of the one-parametric curves in the {h, qs/fi^) —plane 
shown on fig. 4(b). The curves from and S+ and from A_ and S_ lie above and below the 
axis = 0, respectively. The points in which the curves from different families cross each other 
correspond to integer values of both A'^i and N2 and their coordinates define the quantized values 
of h and ^3. Trying to compare the expressions (5.5), (5.11), (5.13) and (5.14) for ^3 with the nu- 
merical results we should keep in mind that we calculated only first few terms of the asymptotic 
expansion of the functions 53 = qs{h; Ni) for small values of integer A^i = 0{hP). Nevertheless, 
a careful examination of fig. 4(b) shows that for all quantized values of ^3 except of those close 
to the degenerate vahie ^3 = agreement is very good. 

Let us summarize the properties of the different families of the curves on fig. 4(b). For fixed 
A^i the function q^{h,Ni) from A+ increases as the conformal weight h grows and for large h it 
approaches the maximal value. For different A^i the following hierarchy holds 



x/27 



> q3{h;0)> q3{h;l)> ...>0. (5.17) 



As conformal weight h decreases, the nonleading terms in the asymptotic expansion of qsih; Ni) 
become important. Using the first 8 terms in the expansion of 93 we find that the function 
q3{h; Ni) vanishes as h approaches the value h = 2Ni + 3 (see fig. 4(b)) 

qsih; Ni)^0, for /i = 2A^i + 3 = Z+ . (5.18) 

We recall, however, that the quantum number qs — has a special meaning as corresponding 
to the degenerate solutions of the Baxter equation (2.11). Notice that equality in (5.18) is 

approximate because the integer Ni takes the vahie. A'^i = 0{h), beyond the approximation 
(5.2), A'^i = 0{h^), under which the corrections to quantized ^3 have been found in (5.14). 
Therefore, using our expressions for q^ih] Ni) we are not allowed to approach 53 = very closely. 
That is the reason why all curves from and should be terminated at the vicinity of the 
axis g3 = 0. To describe the region near ^3 = we have to improve the asymptotic approximation 
for g3 = q3{h;Ni). In contrast with A+, the curves q^ = q3{h;h — N2) from _B+ are decreasing 
functions of h for fixed A^2- They take maximal value q3{N2; 0) at h = N2 and for N2 = 3, 4, ... 
all these points belong to the curve q3{h;0) from A^. As h grows the function q3{h;h — N2) 
decreases and it crosses subsequently another curves from A+ with Ni = 1, 2, N2 — 3. At 
h — 2N2 — 3 = 2A^i + 3 it takes a zero value in accordance with (5.18). 

Although we cannot approach small values of the quantized ^3 within the approximation (5.2), 
a careful examination of fig. 4(b) allows us to speculate about possible behaviour of the function 
Q3(/i; Ni) below the axis ^3 = 0. The comparison of figs. 4(a) and (b) suggests that close to the 
axis ^3 = the different curves from A^ can be considered as a continuation of the curves from 
B_ to the upper half-plane qa > 0. One of such possible curves is shown on fig. 4(a). Together 
with (5.15) and (5.16) this assumption corresponds to the following property of the function 

qsih; Ni) ^ -qs{h; h - 3 - Ni) , < < h - 3 . (5.19) 

This relation is consistent with (5.18). It allows us now to identify the curve on fig. 4(a) as 
Q3 = Q3{h;6). The value of the integer, Ni = 6, can be fixed using (5.18) from the position of 
zero g3 = 0. 
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5.3. Quantization conditions for higher n 

Let us generalize the analysis of the quantization conditions (5.1) to the higher n Reggeon states. 
As we have shown in Sect. 5.1, the possible values of quantized ^3, g„ are parameterized by 
integers A^i, A^„_i. According to the definition (5.4), their values satisfy either Nj — 0{hP), 
or Nj — 0{h). Let us solve the quantization conditions (5.1) in the special case 

ATi , , ... , Ar„_2 = 0(/i°) , Ar„_i = 0{h) , (5.20) 

when almost all roots of the solution of the Baxter equation belong to a single interval. As we 
will see in a moment, the advantage of (5.20) is that for small A"i, Nn-2 one is able to find 
analytical expressions for quantized q^, This does not mean however that the quantization 

conditions (5.1) can not be solved for the values of integers different from (5.20). In the latter 
case the equations for quantized q^, qn become more complicated and their solution is more 
involved. 

Using (5.3) we represent the quantization conditions (5.2) in the following form 
f fix f fix 

where the functions $q and were defined in (3.8) and (4.1). As was explained in Sect. 5.1, for 
each j = 1 , n — 2 the solution of the first condition (5.21) is that the quantum numbers ^3°'', 
qn^ take their critical values, that is they belong to the critical n — 2 dimensional hypersurface 
(3.18). The system of n — 2 equations, (5.21), fixes their position on S„ "almost" uniquely. 
Namely, the solution of (5.21) defines some points on E„, which we identify as quantized values 
of ^3°^, qn^ . For n — 3 we found ^'3°' = ±l/-\/27, but for n — 4 these points can be identified 
as "corners" of the generalized triangle on fig. 5: 

(^^=0, Q'r = ^) and (^^=4' ^^"' = -^)- ^^-22) 

To solve the second equation in (5.21) and find the next corrections, ^3^^ qn\ we have to 
define the limit in which ^3°', qn^ approach their quantized values. Similar to the situation 
for n = 3, we introduce small parameters ^i, £n_2) which measure the size of the infinitesimal 
intervals, Sj = o^j — o^j-x- None of these intervals contains the origin a = and their boundaries 
satisfy the relations Vo(cr2j) = Vo(cr2j_i) = or Vo(cr2j) = Vo(cr2j_i) = 4. Let us concentrate on 
the latter case. Then, the potential Vo(A) has a local maximum on the interval [a2j,a2j-i] at 
a — cr*{ej) and 

<^2,,2,-l = <T*(£,)±|+0(£5). 

The expansion of Vo(A) around A = a* looks like 

Vo{X) = Voia*iej)) + ^K)"(a;(£,))(A - a*)' + ... . 

Substituting A = <72j~i + £jX and taking into account that Vo(A) = 4 for x = and x = 1, we 
obtain the expansion of the potential in the limit Sj — > as 

Vo{a2,-i + e,x) - Voia*) = 'fK'{a*) x{x - 1) + 0{ej) . (5.23) 
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Here, a* = cr*(0) is the position of the maximum of Vo(A) in the hmit Sj = a2j — o'2j-i 0. It 
can be found as a solution of the equations 



VniaV. 







and V^{a*)<0. 



(5.24) 



Repeating analysis for Vo{a2j) — Vo{a2j-i) — one can show that the same expansion (5.23) 
holds provided that we define (7| as a local minimum of the potential, 



VniaT. 











and 



V^'{a*)>0. 



(5.25) 



One can check that for n = 3 the relation (5.23) coincides with (5.7). Substitution of the 
expansion (5.23) into (3.8) and (4.1) yields 



^o{a2j-i + ejz) 
^[{cr2j-i + ejz) 



2 
1 



2\V^'ia*)\Jziz-l) + 0{e 



1 /I 



4\z 



+ 



^z{z-l) V2|K)"« 



where the potential Vi{X) was defined in (2.29) and j — 1, n — 2 numerates the intervals inside 
S. The evaluation of the contour integral of $'i(A) in (5.21) is similar to that in (5.10). It is 
given by the residue of ^[{a2j-i + Sjz) at 2; — > oo and EjZ = fixed. Finally, the solution of the 
second equation (5.21) has the form 



2\VJ/(a*)\ Ni + 



J = l,...,n-2 



Using the definition (2.29) of Vi{X) we obtain from this relation the system of {n — 2) equations 
for quantized ^3^^ qn^ 



f + + (If -1 - Vw«^(iv, + i) . 



(5.26) 



Here, a* is defined as a solution of (5.24) or (5.25). As an example, we consider the system of 
equation (5.26) for the n = 4 Reggeon compound states. Using two different sets (5.22) of the 
quantum numbers (^3°', ^4°') we find the solution of (5.26) in the following form 



-^{N,-N2), 



?4 



1 ^2 
8 ~ ~8~ 



{Ni + N2 + I) 



for the first set in (5.22) and 

g« = ±\{n, - 2V2N2 -I-V2), qi''-^- ^{n, - V2N2) 



for the second set in (5.22). Here, < Ni , N2 < h — 4 are integers parameterizing all possible 
solutions for and q^. 
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6. Asymptotic expansions for the energy of the Reggeon 
states 



In previous sections we developed the method for evaluation of the energy and conserved charges 
of the n Reggeon states in the quasiclassical approximation, h ^ oo. The resulting expressions 
have a form of the asymptotic expansion in powers of 1/h. Although they were found for large 
h, it is of most importance to calculate their values for the conformal weight h = 1/2, which 
corresponds to the maximal energy of the Reggeon state, (2.13), and which defines the intercept 
of the QCD Pomerons and Odderons (1.6). The natural question appears whether it is meaningful 
to substitute h = 1/2 into the large h asymptotic expansion and estimate the approximate value 
of the energy. To test this idea we start with n = 2 Reggeon states and apply (4.16) and (4.23) 
to find the first few terms of the expansion of e2{h) in 1/h 

/,v ,,/,^x2 1 1 1 1 1 691 1 /, ir\ ,„ , , 

s,(h) = -41n(te-.)+-+3^-— +^-g^+33^-5j^+3^+0(ft--) . (6.1) 

One can check that this series coincides with the asymptotic expansion of the exact result (2.19). 
The series (6.1) is a sign-alternating and its coefficients rapidly grow in higher orders. It is 
obvious that for h — 1/2 the series (6.1) becomes divergent while the correct answer for the 
energy £2(1/2) is known to be finite (4.18). On the other hand, as follows from (4.17) and (4.18), 
the sum of the first four terms of the divergent series (6.1) gives for h = 1/2 the result which is 
close to the exact answer. These two properties are remarkable features of the asymptotic series 
[30, 31, 32], which we will explore below to find the energy of the Reggeon states. Our goal will 
be to identify the partial sums which give the best asymptotic approximation to £„. However, 
any asymptotic approximation should be accompanied by bounds for the corresponding error 
terms and in the case of the asymptotic series their determination becomes extremely nontrivial 
[30, 31, 32]. In what follows, instead of giving the rigorous proof we apply a "practical" version 
of the asymptotic approximations [32] based on the well-known examples and briefiy described 
in Appendix C. 



6.1. Euler transformation for n = 2 Reggeon states 

Let us apply the "practical" asymptotic approximation to the asymptotic series (4.16) and (4.23) 
for the energy of the n — 2 Reggeon state and then compare the approximant e2^^{h) with the 
exact result (2.19) for £2{h). To start the analysis, we substitute h — 1/2 into (6.1) and examine 
the first few terms of the series 

£2(1/2) = 0.46 + 4 + 1.33 - 0.53 + 1.02 - 4.27 + 31.03 - 345.58 + 5461.33 + ... . (6.2) 

We find that apart from the first term their absolute value decreases until the n = 4 term is 
reached and then it rapidly increases. This suggests to truncate the series (6.2) and approximate 
£2(1/2) by a partial sum of the first 4 terms. For our purposes, however, it is convenient to exclude 
the 4-th term from the partial sum because after this the approximant has a simple interpretation. 
We notice that the 4-th term in (6.1) and (6.2) has a sign opposite to all preceding terms and 
its inclusion diminishes the partial sum. Therefore, the sum of the first 3 terms corresponds to 
the local maximum, S'max, of the partial sum Sn as a function of a discrete variable n, 

4PP(/i) = 5-^3^ = max^^^_2,...'S'n ■ (6.3) 
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It is this interpretation of the approximant which will be used below. 

Thus defined asymptotic approximation el^^{h) to (6.1) coincides with (4.17) and it is in 
a good agreement with the exact result (4.18). We notice that (4.17) provides us asymptotic 
approximation to S2{h) not only for h = 1/2 but for larger values of the conformal weight h 
as well. The accuracy of the approximation (6.3) and (4.17) is not clear yet. However, the 
fact that the asymptotic series (6.1) and (6.2) are sign-alternating allows us to estimate the 
optimum remainder term [31] (see Appendix C for details). If we assume that the first three 
terms in (6.1) and (6.2) correctly describe the asymptotics of £2{h), then the absolute value of 
the optimum remainder term, Se2{h) = e2{h) — 82^^ (h), is smaller than the last term included 
into the approximant and the first term excluded from the approximant. In application to the 
series (6.2) this property allows us to estimate the remainder as —0.53 < 5e2{l/2) < , which is 
in agreement with the numerical results (4.18). 

Let us apply the Euler transformation [30, 32] (for definition see Appendix C) to find better 
approximation to the maximal energy 62(1/2). To this end, we take the asymptotic series (6.1) 
and reexpand it around the point h + Xq up to 0{{h + Xo)~^^) terms. Although the original 
series (6.1) contains only even powers of 1/h in higher orders, the transformed series has all 
powers of l/(/i + xq). Let us denote the partial sum of the first n terms of the transformed series 
as Sn{h,Xo). According to the "practical" asymptotical approximation (6.3), to find the best 
approximation to e2{h) we have to examine the behaviour of Sn{h, function of n, identify 

the lowest value of n = no at which Sn has a maximum and estimate the best approximation 
to e2{h) as Smax = 'S'„o(^'^o)- To understand the role of the parameter xq in this procedure we 
fix the conformal weight as h — 1/2 and evaluate 5'n(l/2,xo) for n = 1, 15 and two different 
values of xq: xq — 0.25 and xq — 0.5. The results of our calculations are shown on fig. 6. We 




n 



n 



(b) 



Figure 6: The partial sums Sn of the Euler transformed asymptotic series for the energy. (6.1), 
corresponding to = 1/2 and two different values of the parameter of transformation: (a) Xq = 
0.25 and (b) Xq = 0.5. The solid line denotes the exact expression for the energy £2(1/2) = 8 In 2, 
the dotted line indicates the position of the local maximum ^max- 

notice that in both cases the partial sum Sn grows with n until it approaches the local maximum, 
'S'max- Then, the region of stability follows after which Sn starts to oscillate around S'max and the 
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size of the fluctuations rapidly grows with n. It is important to reahze that the position, n, of the 
maximum Sn — -S'max depends on the parameter xq. For xq — 0.25 we have S^ax — 5.489602 at 

n = 5 and for xq = 0.5 it is S'^ax — 5.539171 at n = 10. The dependence of S^ax{h, Xq) on xq for 
/i = 1/2 is shown on fig. 7(a). The cusps correspond to the "critical" values of Xq, at which the 
number of terms contributing to ^'max increases by 1. As parameter xq increases the number 




0.2 0.4 0.6 0.8 2''''i''''i''''i''''t''''-)''''i 

Xo Xo 



(a) (b) 

Figure 7: The dependence of the asymptotic approximation to the energy, S^^^xih = 1/2, a;o), on 
the parameter of the Euler transformation xo for (a) n = 2 and (b) n = 3 Reggeon compound 
state. The dotted hne represents the exact value of the energy e2{h— 1/2) = 81n2. 

of terms contributing to S'max also increases and, hence, the approximation with which one can 
evaluate ^max becomes better. However, for large xq the position of the local maximum of Sn 
is shifted toward larger values of n and in order to apply the Euler transformation we have to 
add higher l//i— terms to the asymptotic expansion of the energy (6.1). Namely, the expression 
(6.1) defines the asymptotic expansion of the energy up to 0{h~^^) terms and under the Euler 
transformation the maximal value of the n corresponding to S'max = Sn cannot exceed n = 15. 
One can show that this condition is still satisfied for xq = 0.75 and the corresponding value of 
the partial sum gives the approximate value of the energy 

4PP(l/2) = 5.545150 

which is remarkably close to the exact answer (4.18). One might expect that the accuracy 
of the asymptotic approximation should be better for larger values of the conformal weight, 
h > 1/2. Indeed, taking the parameter of the Euler transformation as Xq = 0.75 and repeating the 
calculation of S'max for different values of h we obtain the results for 82^^^ (h) and 6e2{h) summarized 
in Table 1 . We conclude that the "practical" version of the asymptotic approximation combined 
with the Euler transformation works perfectly well for the energy of the n — 2 Reggeon states. 
Although we are not able to find the exact value of the energy we have the method which allows 
us to find the approximation, £2^^(/i), and improve its accuracy by varying the value of the 
parameter of the Euler transformation, Xq. This gives us a hope that the same procedure can be 
applied to the evaluation of the energy of the higher Reggeon states. 
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h 


0.5 


1 


1.5 


2 


2.5 


3 


3.5 


4 


app 
^2 


5.5452 


~ 10-5 


-2.4548 


-4.0000 


-5.1215 


-6.0000 


-6.7215 


-7.3333 


862 


~ 10-5 


~ 10-5 


~ 10-^ 


~ 10-^ 


~ 10-8 


~ 10-9 


~ 10-9 


~ 10-9 



Table 1: The asymptotic approximation to e2{h) based on the Euler 
transformation for = 0.75. 



6.2. Analytical continuation 

To apply a similar consideration to the higher Reggeon states we need to know the asymptotic 
expansion of the energy similar to (6.1). In general, it can be found from (4.23), but examining 
the definitions (4.21) and (2.28) we immediately realize an important difference between n = 2 
and n > 3 cases. For higher Reggeon states the energy depends on the quantum numbers 
^3, g„. Their values are quantized and depend on the conformal weight h as well as on the 
set of integer numbers defined in (5.4). The latter dependence was studied in Sect. 5 and the 
expressions for the quantum numbers were obtained in the form of asymptotic series in 1/h. To 
find the asymptotic expansion of the energy of the higher Reggeon states we have to supplement 
the relation (4.23) by analogous expressions (5.4) for the quantum numbers 93, g„. Moreover, 
as we have shown in Sect. 4, the quasiclassical expansion of the energy in powers of I//1 becomes 
divergent. The lowest order terms in the expansion of in powers of 1/h may depend, in general, 
on an arbitrary higher order terms in the expansion of the quantum numbers ^3, g„ and this 
makes the definition of the asymptotic expansion of the energy very problematic. 

Luckily enough, a close examination of the expression (4.23) shows that the 0{h^^)—tcnQ in 
the expansion of the energy e„ depends on the finite number of nonleading terms in the expansion 
(2.27) of the quantum numbers: qn\ Qn-i , q^~^^"~^^\ Therefore, having an expression for the 
quantum numbers up to 0{h~^) terms we will able to derive from (4.23) the expansion of the 
energy with the same accuracy. To find the asymptotic scries for ^3, g„ we have to solve the 
quantization conditions (5.1) and (5.2) as was described in Sect. 5. Their solution is well-defined 
for integer values of the conformal weight h > n and it was parameterized in (5.4) by the set of 
integer numbers A^i, Nn-i- As a result, the energy of the n Reggeon states also depends on 
the same numbers 

En = Snih, gs, Qn) = £n(^; Ni, Nn-2) ■ (6.4) 

This implies, in particular, that similar to the distribution of the quantized ^3, g„, all possible 
values of the energy belong to the families of curves in the {h, £„)— plane labelled by integers Ni, 

Let us find the explicit form of the function £3 = 63{h; Ni) for the n = 3 Reggeon states. 
Substituting the asymptotic expansion (2.27) and (5.14) of the quantized ^3 = qs{h; Ni) into 
(4.23) and (4.21) we find the first 8 terms of the expansion of £3 in powers oi 1/h as 

ss(h; N,) = -6 In j + ^i^^ + ... + -l^^ + 0{h-') , (6.5) 

where the coefficients are given by 

= 6A^i + 6 
ei^ = 5iV,^ + 8iVi + H 
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9 3 ^ 9" 9 

1553 ,- 2336,. 5359,. 190883 

si*' = iV-f + -72Ni 

^ 162 ^ 81 ^ 162 ^ 4860 

.5, 18538 ,,, 13928,,. 32402,. 96188 ,,0 1735909,, 135976 

4'^ = Nf H Nf Nf N? A^i 

^ 1215 ^ 243 ^ 243 ^ 243 ^ 3645 729 

55345 ,. 83152 ,. 625865 ,. 3545360,,, 13670593 ,,0 

e^^' = NY H A^f Nf Nf N, 

^ 2187 ^ 729 ^ 1458 ^ 2187 ^ 4374 ^ 

5748704,, 80780333 



e 



2187 91854 

1983742 ,,7 1489976, ,fi 24244915,,, 111180580,,. 902984221 ,,. 
- n; H iVf N^ iVf Nf 



45927 ' 6561 ' 19683 ' 19683 ' 59049 
398272432 . 2037528641 259171832 



AT^ AT 

19683 ^ 137781 ^ 59049 

The expression (6.5) should be compared with the numerical results for the energy of the n = 3 
Reggeon states presented on fig. 8(a) (see also fig. 3 in [22]). Let us perform the comparison 
following the same procedure as we have used in Sect. 5.2 for quantized qs — qz{h] Ni). Similar to 
(5.15) and (5.16) the function (6.5) defines two different famihes of the curves in the {h, £3)— plane 
shown on fig. 8(b), 

A: £3 = £3(/i; A^i) , /i>2A^i + 3; B : e^^ e^{h]h-Ni-3) , 2N-i^ + 3>h> Ni + 3 , 

(6.6) 

with A^i > being an integer. Notice that due to the invariancc (2.12) under replacement 
— > — ^3 the energy e^{h\Ni) takes the same value for the curves from and A_. The 
quantized values of the energy correspond to the points on fig. 8(b), at which the curves from 
A and B cross each other. There is one-to-one correspondence between different curves on 
figs. 4(b) and 8(b), such that each function ^3 = q^lh; Ni) is mapped into the corresponding 
function £3 = e3{h; Ni). As an example, the solid line on fig. 4(a) corresponding to ^3 = qsih; 6) 
induces the curve £3 = £3(/i;6) shown on fig. 8(a). In the same manner, the axis ^3 = on 
fig. 4(a) is mapped into the dashed line on fig. 8(a) describing the degenerate values of the 
energy es{h, qs — 0) = S2{h). Similar to the distribution of the quantized ^3 on fig. 4, we observe 
a good agreement on fig. 8 between numerical results for the energy and analytical expression (6.5) 
everywhere except of the points close to the "degenerate" curve £3 = e^lh, qs = 0). Nevertheless, 
as one can see from fig. 8(a), the function e3{h,q3) is a smooth function of h at the vicinity of 
qs — O and according to (5.19) it satisfies the following relation 

e3{h;m)^S3{h;h-Ni-3), (6.7) 

which allows us to consider (6.6) as definition of the unique function £3 = £3(/i; A^i) for A^^i + 3 < 
h < 00. Solving the Baxter equation, we are interesting to find the maximal value of the energy 
of the n Reggeon states, (1.6) and (2.8). As was shown in [22] and as it can be easily seen from 
fig. 8(a), for positive integer conformal weight, h > n, the energy Cn cannot exceed the upper 
bound 

en{h,qs,...,qn)<-^^{h)-jjj{l)), h>n (6.8) 

given by the energy of n = 2 Reggeon state with the conformal weight h. As h decreases the 
maximal value of the energy max^^ ^^ £„(/?., ^3, g„) grows until h approaches the value h = n, 
beyond which the Baxter equation does not have polynomial solutions. On the other hand, we 
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h h 

(a) (b) 

Figure 8: Holomorphic energy of the n = 3 Reggeon compound states: (a) the results of the 
numerical solution; (b) two different famihes of curves defined in (6.6). The solid hne on (a) 
corresponds to £3 = S3{h, qs) with the function — qz{h] 6) shown on fig. 4(a). The dashed line 
on (a) represents the energy of the n — 2 Reggeon states (2.19) and the dotted line on (a) is 
described by the function £^{h] 0). 



expect from (2.13) that the absolute maximum of energy 8^ is achieved for h = 1/2 and not for 
h — n. This means that in order to approach physically interesting values of the energy we have 
to penetrate through the h — n barrier to the region of smaller values of the conformal weight 
h. In what follows we restrict consideration to the n — ?> Reggeon states while generalization to 
the higher Reggeon states is straightforward. 

To perform analytical continuation of the energy e^ih, q^) for /i < 3 we use fig. 8(a) to observe 
the following properties of the function e^{h; Ni). Let us consider the flow of the quantized values 
of the energy along the curve £3 = e^{h] Ni) for fixed Ni. The energy increases along es{h; Ni) 
as conformal weight decreases from h = 00 to h = 2Ni + 3. At the point h = 2Ni + 3 it takes the 
degenerate value e3{h; Ni) = e2{h) corresponding to qsi^h; Ni) = in (5.18). For smaller values 
of the conformal weight, h < 2Ni + 3, the curve S3 = e3{h;Ni) is "reflected" from the upper 
limit £3 = e2{h) and ior h — Ni + 3 it approaches the end-point, £3 = Ss{h; Ni) — es{h; 0), and 
is terminated. 

We notice from fig. 8 that any two curves e-i{h; Ni) and e-i{h; N[), corresponding to different 
values of integers A^i and N[, cross each other at only one point. Its position can be easily found 
from (6.7) as 

esih; Ni) = e3ih; N[) , for h = Ni + N[ + 3. 

It is clear, that if the value of one of the functions is bigger for h < iVi + iV( + 3, then for 
h > Ni + N[ + 3 the relation between them becomes opposite. As a result, the functions satisfy 
the following hierarchy 

-4{ip{h) - > £3(/i; A^i) > £3(/i; N[) , for h>Ni + N[ + 3, Ni > N[ 
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and 

-4{ij{h)-ij{l))>e3{h;N[)>es{h;N,), for h < + N[ + 3 , N,>N[. 

These inequalities allow us to identify the maximal value of the energy corresponding to the 
integer conformal weight h > 3. We find that for the values of the conformal weight inside the 
interval 2k + 2 < h < 2k + A, with k positive integer, among all the functions £3 = £3(/i; A^i) the 



h-2 



takes a maximal value 



one with Ni — k 

max^^esih; Ni) = e3[h; [^Jj . (6.9) 

Therefore, as h decreases from h = 00 to h = 4 the maximum of the energy "jumps" from one 
curve S3{h; Ni) to the next one, A^i —>■ Ni — 1, and finally for /i < 4 it follows the function £3(/i; 0). 

For h > 3 the quantized values of the energy are confined to the region on fig. 8 which is 
restricted from above by the function £3^/1; j and by the function £3(h; 0) from below, 

mm^^e3{h;N,)=es{h;0). (6.10) 

As conformal weight h decreases, the upper and lower bounds for the energy, (6.9) and (6.10), 
respectively, move toward each other leaving no phase space for quantized values of energy £3. 
At h < A both boundaries coincide and the energy is given by 

£3 = £3(/i;0), for /i<4. (6.11) 

For h — 3 it takes the value £3 = £3(8; 0) = —6 corresponding to the n — 3 degenerate Reggeon 
state with q3{3; 0) = —0. It is of the most importance now to understand the behaviour of the 
function £3(^1; 0) beyond the barrier, h < 3, where the Baxter equation does not have polynomial 
solutions. We conclude from (6.9) and (6.10) that once the conformal weight approaches the 
barrier h = 3, among all functions S3{h;Ni) describing quantized values of energy only one, 
(6.11), with A^i = survives. Therefore, performing analytical continuation to /i < 3 we may 
restrict ourselves to the function es{h; 0). 

6.3. Euler transformation for n = 3 

To define £3{h; 0) for /i < 3 we use its asymptotic expansion in powers oil/h and apply the Euler 
transformation to evaluate the best asymptotic approximation. The asymptotic expansion for 
g3(/i;0) can be found from (2.27) and (5.14) up to 0{h~^^) terms as^ 

\/3, H / 3 4 8 688 7360 

261472 10272032 447738400 7244437856 .^/,-io\\ 
~6561/i6 ~ 59049 /i^ ~ 531441 /i^ ~ 1594323 /i^ ^ ^ H' 
As was explained in Sect. 6.2, substituting (6.12) into (4.23) we can obtain the expansion of the 
energy with the same accuracy 

(he^^\ 6 11 56 190883 135976 

80780333 259171832 506229906961 226978281368 ,^.,_lo^ 



91854 59049 /i^ 21257640 /i^ 1594323 /i^ 



^In comparison with (5.14) we calculated and added to q3{h;0) two additional nonleading terms. 
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Let us compare this expression with an analogous expansion (6.1) of the energy of the n — 2 
Reggeon states. We notice that in both cases the absolute value of the coefficients in front of 
powers of rapidly increases in higher orders revealing the asymptotic nature of the series. 
However, in contrast with (6.1) the series for e3{h;0) is not a sign-alternating. Starting from 
0{h~^) all terms in the series (6.13) are negative and this makes its properties different from 
(6.1). As one can see from fig. 8(b), the expression (6.13) is in a good agreement with the 
numerical results for the energy for large conformal weight h > 5. For small h the absolute value 
of the higher order terms in (6.13) starts to grow and the series (6.13) becomes divergent. In 
particular, substituting h — 1/2 into (6.13) we obtain the series 

£3(1/2; 0) = 3.99 + 12 + 7.33-49.78- 628.42 -5968.77 (6.14) 

-56284.34 - 561804.51 - 6096389.64 - 72891678.83 + ... , 

which should be compared with an analogous series in (6.2). 

To evaluate the energy £3(/i;0) for /i < 5 we have to define the asymptotic approximation 
to the series (6.13). Let us denote the partial sum of the first n terms of the series (6.13) as 
Sn- Then, it follows from (6.13) and (6.14) that the partial sum Sn grows as n increases, at 
n = 3 it takes a maximal value, S'max, and for n > 4 it rapidly decreases. It is clear that S'max 
overestimates the energy while the partial sum Sn for large n underestimates it. This means that 
the remainder term, as a function of a discrete n, first decreases for small n until it approaches 
a minimal negative value corresponding to Sn — -S'max and then it rapidly grows to 00. To our 
practical needs, we identify the asymptotic approximation to the energy as 

ef^{h; 0) = ^max = max„^^ 2 .Sn . (6.15) 

Under this choice, the remainder term Se^ = £3 — £3^^ is always negative and the approximant 
(6.15) defines an upper limit for the energy e^Xh; 0). For the series (6.13) and (6.14) the expression 
for S'max is given by the sum of the first three terms. To increase the number of terms contributing 
to -S'max and, as a consequence, to improve the accuracy of the asymptotic approximation (6.15) 
we apply the Euler transformation to the series (6.13). The behaviour of the partial sum, Sn{h — 
1/2, Xq), as a function of n is shown on figs. 9(a) and (b) for two different values of the parameter 
of transformation, xq = 4 and xq = 7.85, respectively. Changing xq we find from fig. 9 that as 
Xq increases, the growth of Sn with n becomes less steep and the value of S'max slowly decreases. 
Since Smax defines the upper bound for the energy £3(/i; 0), the latter property implies that the 
accuracy of the asymptotic approximation (6.15) becomes better. The dependence of Smax(^; Xq) 
on Xq for = 1/2 is shown on fig. 7(b). It suggests that approximant Smaxj being a decreasing 
function of Xq, asymptotically approaches the energy e3{h; 0) from above for large enough Xq. 

Increasing the value of xq we have to check that the number of terms contributing to Smax is 
less than 10, that is the number of terms in the expansion (6.13). Otherwise, additional higher 
l//i— terms should be taken into account in the asymptotic expansion (6.13). We found that 
the maximal value of Xq for the series (6.13) which satisfies this condition is close to Xq = 7.85. 
Applying the Euler transformation to the series (6.13) and taking Xq = 7.85 we obtain the 
asymptotic approximation to the energy of the n = 3 Reggeon states, sl^^{h;0), shown on 
fig. 10. 

Let us summarize the properties of the function e^^{h; 0). According to the definition (6.15), 
this function approximates the energy of the n = 3 Reggeon states £3 = e^lh; 0) for an arbitrary 
value of the conformal weight 1/2 < h < 00. Although its accuracy was not defined yet, we know 
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(a) 



n 



(b) 



n 



Figure 9: Partial sums 5'„ of the Euler transformed series for the energy of n = 3 Reggeon states, 
(6.13), corresponding to h — 1/2 and two values of the parameter of transformation: (a) xq — A 
and (b) xq — 7.85. The dotted line indicates the position of the asymptotic approximation Smax- 



in advance that the asymptotic approximation (6.15) provides upper bound for the energy. For 
large values of the conformal weight, h > 5, the expression (6.13) is in a good agreement with 
the numerical results for the energy and the asymptotic approximation e^^^{h; 0) makes it even 
better. For smaller values of h the expansion (6.13) becomes divergent and it is approximated 
by the function e^^^{h; 0). To estimate the accuracy of the approximation we use the numerical 



results for the energy, e3{h; 0), and compare them with the corresponding values of e^^^{h; 0) to 
find the remainder terms &3 = £3 — £3^^. The numerical data are summarized in Table. 2. As 
was expected, the remainder terms are negative. 



h 


3 


4 


5 


6 


7 


8 


9 


10 


app 
^3 


-5.0809 


-7.1050 


-8.6696 


-9.9333 


-10.9880 


-11.8907 


-12.6783 


-13.3763 




-0.9191 


-0.3950 


-0.1637 


-0.0667 


-0.0273 


-0.0114 


-0.0050 


-0.0022 



Table 2: The asymptotic approximation to S3{h;0) based on the Euler transformation for xq — 
7.85. 

Performing analytical continuation of the energy to /i < 3 we use the symmetry (2.12) of the 

Baxter equation under the replacement h ^ 1 — h to restrict the value of the conformal weight 
to the fundamental domain, 1/2 < h < 00. We notice from fig. 10 that e^^^{h; 0) is a decreasing 
function of h. This implies that, in accordance with (2.13), the maximal value of the energy of 
the n — 3 Reggeon state corresponds to /i = 1/2 and it is given by 



max, £ 



app 



ih;0) 



4.826208 . 



(6.16) 



We stress that this expression should be considered as upper bound for the energy £3(1/2; 0). 
The function e'^^^{h; 0) depends on the parameter of the Euler transformation xo- For xq = 7.85 
we find from fig. 10 that el^^{h;0) crosses e2{h) = —4{ilj{h) — V'(l)) two points, = 4.38 
and h- = 0.56. For fixed h, the value of £3^^(/i;0) decreases as Xq increases (see fig. 7(b)) and 
the distance between and /i_ becomes smaller. One might expect that the distance /i+ — /i_ 
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2 4 6 8 10 



h 

(a) (b) 

Figure 10: Holomorphic energy of the n = 3 Reggeon compound states: (a) Dots represent 
numerical results for the energy £3 — esi^h; 0). Dotted line corresponds to the asymptotic series 
(6.13). Thick line on (a) and (b), £3 — s'^^'^{h; 0), was found after the Euler transformation of the 
series (6.13) with the parameter xq — 7.85. Thin line on (b) is the energy of the n — 2 Reggeon 
states, S2{h), defined in (2.19). 



controls the accuracy of the asymptotic approximation and in the limit h^ — h^ — > the function 
el^'^{h; 0) coincides with the exact expression for the energy, £3(^1; 0). We recall that for integer 
Ni > 1 the function es{h; Ni) crosses e2{h) at only one point h = 2Ni + 3. If one assumes that 
the same property holds for A^i = 0, then h+ — h- — 3 and es{h;0) < — 4('0(/i) — '(/'(I))- 
h = 1/2 this leads to the relation 

£3(l/2;0) < 81n2, (6.17) 

which is in agreement with (6.16). Finally, taking into account (6.16) we find from (1.6) and 
(2.8) the asymptotic approximation of the intercept of the Odderon as 

a N 

^app ^i^^ 2.413104 . (6.18) 

Odderon ^ ^ ' 

This expression is smaller than the intercept of the BFKL Pomeron, (1.1), and it approximates 
Q^oddoron fro™ above. The estimate (6.18) is also in agreement with the lower bound for the 
Odderon intercept proposed in [11]. 

One may try to repeat the analysis and apply the Euler transformation to the asymptotic 
series (6.12) to find the approximation to the quantum number 53(1/2; 0) corresponding to the 
maximal energy of the n = 3 Reggeon states. We find however that since all terms in the power 
series (6.12) except of the 0{h~'^) term have the same sign, the Euler transformation does not 
improve its convergency properties and we are not able to construct the convergent sequence of 
the partial sums similar to (6.15) for the energy. Moreover, as we will show in the next section, 
the coefficients in (6.12) grow as factorials and the asymptotic series for g3(/i;0) is not Borel 
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summable. Unfortunately, there is no regular way for approximating non Borel summable series 
like (6.12) and in order to calculate ^3(1/2; 0) one has to develop another summation method. 



6.4. Borel summability 

As was found in the previous sections, the asymptotic expansion of the energy of the n— Reggeon 
states has the following general form 

Snih, {q}) = -21n(/i"g„e"^-) + ^ ^ (6.19) 

k>i 

where q2 = Ij ^3 = ■■■ and f\ are rapidly growing coefficients. It is important now to 

understand whether the asymptotic series (6.19) defines uniquely the energy. To this end we 
perform the Borel transformation of the power series entering into (6.19) 

Snih, {q}) = -2 ln(/i"g„e"^K) + / dt B,^{t) e~^' . (6.20) 

Jo 

Here, the function B^^it) is defined as 

BeSt) = E ^r^'^ = E ^M.} e-^* - ^ (6.21) 

fe>0 fe>0 

where in the last relation we used (A. 9) and (A. 11) and the additional term, —2n/t, was added 
to provide regularity of B^^ as t — > 0. The properties of the function B^^{t) depend on the large 
order behaviour of the coefficients fk- Let us introduce the following ansatz [33] 

h^ck\a-^k^{\ + 0{l/k)) as A; ^ oo , (6.22) 

which being substituted into (6.21) leads to the singularity of the function B^^it) ai t = a 

(6.23) 

To verify (6.22) one has to compare fk with their numerical values in the asymptotic expansions 
(6.1) and (6.13) and extract the values of the parameters a, b and c. 

For the n = 2 Rcggcon states one uses the last three terms in the expansion (6.1) to find the 
following values: a = ±6.2746 i, b = —1.0201 and c = —8.2657, where the sign of a can not fixed 
since the expansion (6.1) contains only even powers of h in higher orders. These values are very 
close to the exact expressions 

a = ±27ri, b = -1 , c = -8 , (6.24) 

which follow from the well-known asymptotic expansion of the ^—function in (2.19). Moreover, 
the function B^^{t) can be evaluated explicitly using (6.21) and (A. 10) as 

4 4 
B,Jt) = 

We check that, in accordance with (6.23) and (6.24), the nearest to the origin singularity of this 
function is located at t = ±27ri. At the same time, the function B^^{t) does not have singularities 
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on the integration path in (6.20) and as a consequence the asymptotic series (6.20) for the energy 
of the n — 2 Reggeon states is Borel summable and it can be approximated using the well-known 
methods [33]. 

Let us check the large order behaviour of the energy of the n — 3 Reggeon states (6.13) using 
the ansatz (6.22). We find that for the set of parameters 

a = 1.8098, 6 = 1.8195, c = -1.5000 (6.25) 

the relation (6.22) reproduces the coefficients /4, /g with 8% accuracy. This result indicates 
that the function Bs^{t) has singularity (6.23) at the point t — 1.8098 which belongs to the 
integration path in (6.20). Therefore the Borel transformation does not exist and the asymptotic 
series for the energy of the n = 3 Reggeon states is not Borel summable. Different prescriptions 
for integrating the singularity of B^,^{t) lead to the results for the energy which differ in ~ 
exp(— /la) terms. Thus, the asymptotic series (6.13) does not define the energy uniquely but rather 
represents one member of an infinite class of functions having the same asymptotic expansion. 
To avoid ambiguity and fix a particular member of the class we have to provide additional 
information about analytical properties of the energy. This requires solution of Baxter equation 
beyond the large h approximation. 

We notice that the asymptotic expansion of q^{h] 0) defined in (6.12) exhibits the large order 
behaviour similar to (6.13). Indeed, rewriting (6.12) as 

^3(/^;O) = ^/.-^(|l + E0j (6.26) 

we find that the coefficients {k = 4,..., 9) satisfy the asymptotic behaviour (6.22) with 2% 
accuracy provided that 

a = 1.8174, 6 = 0.8019, c = -0.4650 . (6.27) 

As a consequence, the asymptotic series for q3{h; 0) suffers from the same problems as that for 
S3{h; 0). We notice that the expressions (6.25) and (6.27) suggest the following relation between 
asymptotic expansions for the n — 3 Reggeon states: guk/ fk ~ 0.28. 



7. Conclusions 

The spectrum of the color singlet n Reggeon compound states - perturbative Pomerons and 
Odderons, is expressed by means of the Bethe Ansatz in terms of the solution of the Baxter 
equation for the XXX Heisenberg magnet. In this paper wc developed the method which allows 
us to find the solution of the Baxter equation for an arbitrary number of the Reggeized gluons. 

The method is based on the observation that in the limit of large integer conformal weight h of 
the Reggeon states, the Baxter equation takes a form of a discrete one-dimensional Schrodinger 
equation for a particle in the external potential ^(A), which depends on the conserved charges 
^3, g„ of the Reggeon states and is singular at the origin A ^ 0. The inverse conformal 
weight, plays a role of the Planck constant and this fact allows us to apply the well-known 
quasiclassical expansion and obtain the solution of the Baxter equation, $(A), as an asymptotic 
series in 1/h. We found that the derivative of the solution, $'(A), has a discontinuity across a 
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finite disconnected interval <S on the real axis in the complex A— plane, whose position and the 
number of connected intervals inside S depend on the potential or equivalently on the value of 

the charges ^3, g„. The quantization conditions for the charges can be reformulated as certain 
constraints on S. Namely, for the n Reggeon compound states the interval of nonanaliticity 
S should consist of n — 1 connected intervals. The corresponding quantization conditions have 
solutions for charges which can be parameterized by the set of integer numbers A^^i, Nn-i. 
We have shown that the value of quantized qs, Qn can be obtained as a series in and the 
expressions for the lowest order coefficients were found. We observed that for n = 3 and n = 4 
Reggeon states they are in a complete agreement with the results of the numerical solutions of 
the Baxter equation. 

The energy of the n Reggeon states depends on the conformal weight h and quantized charges 
qs, Qn- Using the quasiclassical solution of the Baxter equation we obtained the expansion of 
the energy in powers of 1/h for h>n and then analytically continued the result for small values 
of conformal weight, 1/2 < h < n. For large integer h it agrees with the results of the numerical 
solutions while for small h the asymptotic expansion of the energy becomes divergent and it 
should be replaced by the asymptotic approximation. Using the properties of the asymptotic 
series and applying the Euler transformation we defined the asymptotic approximation of the 
energy for small h and checked its validity in the special case of the n = 2 Reggeon states by 
comparing with the exact expression for the energy. As first nontrivial application, we derived 
the asymptotic approximation of the energy of the n = 3 Reggeon state, perturbative Odderon, 
and estimated its intercept. The generalization of the obtained results to the higher Reggeon 
states is straightforward. 

This work was supported in part by the National Science Foundation under grant PHY 
9309888. 



Appendix A. Analytical properties of the energy 

Let us show that the energy of the n— Reggeon compound states satisfies the conditions of the 
Carlson's theorem and therefore it has the unique analytical continuation in the region (2.14). 

We start with the simplest case of the n = 2 Reggeon states, in which the holomorphic energy 
depends only on the conformal weight h and it is given by (2.19). The function e2{h) has poles 
at the origin h — and at negative integer h 

£2(/i)'~"'^, A; = 0,1, 2,... (A.l) 
and its asymptotics at infinity is 

e^{h)''^°° -4\nh. (A.2) 

Thus, the function e2{h) has the unique continuation from integer to complex h defined in (2.3). 

Let us generalize the relations (A.l) and (A.2) to the n— Reggeon states. Their holomorphic 
energy was obtained in Sect. 4 in the form of the asymptotic series inl/h whereas the analytical 
expression for similar to (2.19) is not available yet. In the large h limit the asymptotic 
behaviour of £„ was found in (3.22) as 

en{h)^^°° -2nlnh, (A.3) 
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and it is in agreement with (A. 2) for n — 2. Trying to satisfy the second condition of the Carlson's 
theorem we find that having only a few first terms, (4.16) and (4.23), of the asymptotic series 
for En we are not able to identify all singularities of in the complex /i— plane. 

There is however another way to study the singularities of the energy of the Reggeon com- 
pound states based on their relation with the anomalous dimensions of composite higher twist 
operators entering into the operator product expansion of the structure function of deep inelastic 
photon-hadron scattering at small values of the Bjorken variable x [5] . Let us consider the struc- 
ture function at small x in the generalized leading logarithmic approximation [22]. We choose for 
simplicity the hadron to be a perturbative onium state built from two heavy quarks and created 
via the decay of the photon with invariant mass m. The structure function F{x, Q"^) is defined 
as an imaginary part of the amplitude of the forward photon-onium scattering with the center- 
of-mass energy s — Q'^{1 — x)/x and photon virtuality —Q^. Let us consider the asymptotics of 
the structure function in the double scaling limit of large Q"^ /m^, and small x ~ j s. 

In the limit ^ m^, the structure function can be expanded in powers of us- 
ing the operator product expansion (OPE). In the leading In Q— approximation, <^ \ and 
a^\vi.{Q^ I rr?) ~ 1, this expansion looks hke 

FM) ^ I dxx-'F{x,Q') = ^E^if^(^j ' (A-4) 

where the A;— th term is associated with the contribution of twist— 2/c operators. Their matrix 
elements have the anomalous dimensions 7^^^^^ = 0{as) which have a perturbative expansion in 
powers of ctg. 

In the small x limit, or equivalently in the limit of large energy s ~ Q^/x ^ Q^, the 
structure function has the Regge behaviour which it is governed by the compound Reggeon 
states propagating in the t— channel 



i^.(Q') = i:"rXn(g^), 

n=2 

where the contribution of the n— Reggeon states can be found using (1.3) for s Q'^/x and t — 



as 



{q} ^ -^n,{g} 

Here the summation is performed over all quantum numbers h, q^, qn corresponding to the n 
Reggeon compound state with the energy En,{q}- The residue factors have been defined in (1.5) 
and in the case of the forward scattering, t = 0, they depend only on the invariant masses of 
scattering particles, and m^. In the generalized leading logarithmic approximation, one may 
calculate the residue factors for perturbative states of virtual photon, 7*(Q^), and onium, h{rn?), 
in the Born approximation and neglect corrections. As a result, f^^^Xy* and /^^^/^ do not have 
anomalous dimension and their scaling dimensions are equal to the sum of the scaling dimensions 
of the n Reggeon state, h + h = 1 + 2iz/, and the scaling dimensions of photon and onium states 

fQ{<l} (r)\ — Ml} f^-l+2iu o{q} /r|\ _ Ml} ^-l+2iv 
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where the dimensionless coefficients depend on the quantum numbers of the Reggeon states and 
scattering particles. Substituting these relations into (A. 5) we obtain 




gstr^n^— ^0 a; - ^ [En{^ + {q}) + - iv- {q}) 

(A.6) 

Here, we extracted the sum over quantized values (2.3) of the conformal weight, that is summation 
over discrete m and integration over continious from the sum over all quantum numbers in 
(A.5). 

Let us consider (A.6) in the limit S> m^, in which one should be able to reproduce l/Q^ 
expansion (A. 4). For 3> one can enclose the integration contour over v into the lower 
half-plane, Imi/ < 0, and calculate the integral over v in (A.6) by taking the residue at the values 
of V which satisfy the relation 

= + W )+^n [-^ {q} j ■ (A.7) 

Solving this equation one can find the values of iu which determine the power of m/Q in the 
l/Q— expansion of the structure function (A. 4), or equivalently define the scaling dimensions of 
the composite operators entering into the OPE. The comparison of (A.6) with (A. 4) requires 
that the solutions of (A.7) should have the following form for Imi/ < 0: 

ii^ = -(^ + ^)+C»(«s), Z = 0, 1,.... (A.8) 

Let us now take into account that in order for the contribution of the pole in i/ to (A.6) to be 
nonvanishing, the residue factors, or equivalently the coefficients C^^^* and C|!^/j, should be 
different from zero. This condition imposes selection rules on the quantum numbers h, q^, g„ 
of the Reggeon states. In particular, for the residue factors to be scalar, the conformal spin of 
the Reggeon state, h — h = m, should be equal to the conformal spin of the photon and onium 
states [5]: m = or m = 2. 

Then, substituting (A.8) into (A.7) and taking the hmit ccg — > we find that the holomorphic 
energy Sn = £n{h; {q}) has simple poles at the origin and at the integer negative values of the 
conformal spin 

enih;{q})'---'^, ^ = 0,1,2,... 
n + k 

Moreover, using the asymptotic behaviour (A. 3) and calculating the discontinuity of the energy 
at the negative h one can write the dispersion relation for the function £„(/i; {q}) in the complex 
/i— plane which leads to 

with C some infinite /i— independent subtraction constant. One can check that this relation holds 
for the energy of the n — 2 Reggeon states (2.19) provided that 

Afc = 4 (A.IO) 

and the constant C can be found using the condition £2(1) — 0. The asymptotics of the energy 
(A. 3) implies the following relation: 

74fe,{5}~2n, as /c ^ 00. (-^-H) 
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Thus, the consistency of the small x asymptotics and the large expansion of the structure 
function allows us to identify singularities of the holomorphic energy of the n— Reggeon states 
and, as a consequence, uniquely define £„ in the half- plane Re /i > 1/2 by its values at the integer 
positive h. 



Appendix B. Conformal operators 

The n Reggeon compound states corresponding to the polynomial solutions of the Baxter equation 
have a simple interpretation in terms of the so-called conformal operators [34]-[39]. Let us 
consider the holomorphic wave function of the Reggeon states defined in (1.9). As was explained 
in Sect. 1, the same wave function diagonalizes the hamiltonian of the XXX Heisenberg magnet 
of spin s = 0. The explicit expression for ^n{{zk}') Zq) in terms of the solution of the Baxter 
equation was found in [15, 22] and it is based on the correspondence between eigenstates of the 
XXX Heisenberg magnets of spins s = and s = — 1 

Vn{{Zk}-, Zq) = {ZuZ23--Znl) V^n^~^\{Zk}; ^o) (B.l) 

where are the holomorphic coordinates of the Reggeons and Zq is the holomorophic coordinate 
of the center-of-mass of the Reggeon compound state. Invariance of the Reggeon hamiltonian 
under the conformal transformations (2.2) leads to the following constraints 



(T. ^kdk + zodo + h + r?j <^(f=-^) = , (B.2) 

4dk + 'izk + zld^ + 2hz^ (/^(f =-1) = . 

The expression for </7n({-2^fe}; -^^o) proposed in [15, 22] satisfies these relations. However, instead of 
using the Bethe Ansatz solution for (pn{{zk}\ZQ) we would like to interpret (B.2) as conformal 
Ward identities for {n + 1)— point correlation functions in some two-dimensional conformal field 
theory [5, 14] 

^^r~^\{zkY zo) = {(l){z,)..4{z,,)0h{z^)) . (B.3) 

Here, the field 4>{zk) describes the Reggeon with the holomorphic coordinate z^ and the opera- 
tor Oh{zQ) extrapolates the compound n Reggeon state. To satisfy (B.2), the fields (t){zk) and 
Ohizo) should be quasiprimary operators with conformal weights 1 and /i, respectively [24]. The 
conformal Ward identities fix two- and three-point correlation functions up to a normalization 
constant as 

(0(^1)0(^2)) = const. X z'^2 , {(l>{zi)(l){z2)Oh{zo)) = const, x z^q^ z^q^ z'^^'^ . 

Substituting the last relation into (B.3) and (B.l) one can find the holomorphic wave function 
of the n — 2 Reggeon compound state [3] . 

For integer h > n one can construct Oh{zo) as a composite operator built from n Reggeon 
fields (j) and their derivatives. The corresponding construction has been developed in QCD many 
years ago [35, 36, 37, 39] within the framework of the so-called "coUinear" conformal group and it 
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can be trivially generalized to the "transverse" conformal group which is a symmetry group of the 
Reggeon states. The properties of the resulting conformal operators can be summarized as follows. 
For positive integer conformal weight h >n, one can define the basis of the conformal operators 
and then represent all possible quasiprimary operators Oh{zo) as their linear combination. For 
fixed number of the Reggeons, n, the conformal basis consists of the operators Oi-^^,,,i^_-^{zo) labelled 
by integers h, In-i satisfying the following conditions 

n-1 

h,...,ln-i>0, J2^k^h-n. 

k=l 

The explicit form of the operators Oii,...,/„_i(2;o) for = is [37, 39] 

Oi,..,^_M = a'^-"P,,..„_, f I, %iVi(o)-</'n-i(o)0„(o) , (B.4) 



where d = di + ... + dn, the derivative dk acts on the free field (p^ = (p and -Pii,...,i„_i(ai, On-i) 
is the so-called conformal polynomial in ai, q;„_i. For different values of the 
functions -Pii,...,/„_i(Q;i, «n-i) form the system of biorthogonal polynomials on the simplex 
«! + ... + a„_i + a„ = 1 with the weight function and the measure 

w{ai,...,an-i,an) = {2n-l)\ai...an-ian, J Va ^ dai...dan6(^ at - ij ■ 

In the special case of the n — 2 Reggeon state, thus defined conformal polynomial coincides 
with the Gegenbauer polynomial [34, 35, 36], Pi{a) = Cf^^{2a — 1), and the conformal operator 
interpolating the BFKL Pomeron for integer conformal weight h >2 has a form 

For n = 3 Reggeon states, the conformal operators (B.4) are parameterized by two integers, 
h, h > and h + h — h — 3. The conformal polynomials, Pii,i2, are given by the Appell 
polynomials [36] and the basis of conformal operators looks like 

Oi„iM = d'-' Ji,,i, (^5, 2, 2; |, | j 01 (0)02(0)03(0) (B.5) 

with d = 81+82+03. We recall that an arbitrary conformal operator Oh{0) is a linear combination 
of the basis operators Oi-^ i.^{0). However, for the conformal operator 0^(0) to interpolate the 
n = 3 Reggeon state, the correlation function (B.3) has to satisfy the additional condition of 
being the eigenstate of the conserved charge qs defined in (2.1). Substituting the operators (B.5) 
into (B.3) and (B.l) one can define the basis of the wave functions v^^j./j, in which the operator ^3 
can be represented as a finite-dimensional matrix {q3)i-^,i2;i[,i'^- The eigenvalues of this matrix give 
the quantized values of the charge 93 and the corresponding eigenstates define the decomposition 
of the interpolating operator over the basis of conformal operators (B.5). It is clear that quantized 
values of qs can be parameterized in the same form as the diagonal elements of the matrix 

^3 = ^3(^1,^2) , ^1, ^2 > h + k^ h-3. 
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Comparing this relation with (5.4) we find that they are identical provided that we identify 
integers as A^i = Zi and N2 — I2 + 3. It can be easily seen that the same correspondence holds 
for the higher Reggeon states, A^i = Zi, Nn-2 — ln-2 and A^n-i — In-i + ^■ 



Appendix C. "Practical" asymptotical approximation 

Suppose f{z) is a function of z having the asymptotic expansion 

f{z) = ao + ^ + ^ + --- . (C.l) 

z z^ 

Let us assume that for a given z the absolute values of the first n successive terms of the series 
decrease and then they increase starting from (n + 1)— th term. Then, the asymptotic expansion 
is expecting to give the best approximation to f{z) when it is truncated at about n = n{z) terms 
[32]. The corresponding remainder term 

Rr^iz) ^ fiz) - E 



k=0 ^ 

is called the optimum remainder term. The accuracy of the approximation is controlled by the 
value of Rn{z). 

Euler transformation. 

There is a simple method, the so-called Euler's transformation, which allows us to increase the 
accuracy of the asymptotic approximation [30, 32]. Let us perform the identical transformation 
of the series (C.l) by expanding its each term in a series of inverse powers oi z + xq as 

/(^) = E 3TT = E (.^Uk+1 ' = E J 4 

with xq being a parameter of the transformation. It turns out that under proper choice of the 
parameter xq the sum of the first n' — n'{z, xq) terms of a new series converges more rapidly and 
the corresponding optimum remainder term becomes smaller then for the original series. 

Estimation of the remainder terms. 

Suppose that Sn — J2k=Q is the partial sum of sign-alternating scries = IZfcLo « and Rn is 
the corresponding remainder term. Then, Sn overestimate the Soo when the last term included 
is positive, a„ > 0, and underestimate it when a„ < 0. This means that the remainder term 
has a sign opposite to a„ and it follows from the relation i?„ — Rn+i = On+i that 

\Rn\ < \(^n \ -I \Rn\ < |c^n+l| ■ 

Thus, the absolute value of the optimum remainder term is smaller than the last term included 
into the approximant and the first term excluded from the approximant [31]. 
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